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Abstract 

Morphogenesis, as it is understood in a wide sense by Rene Thorn [l], is considered for various types of chaos. 
That is, those, obtained by period-doubhng cascade, Devaney's and Li-Yorke chaos. Moreover, in discussion 
form we consider inheritance of intcrmittcncy, the double-scroh Chua's attractor and quasiperiodical motions 
as a possible skeleton of a chaotic attractor. To make our introduction of the paper more clear, we have to 
say that one may consider other various accompanying concepts of chaos such that a structure of the chaotic 
attractor, its fractal dimension, form of the bifurcation diagram, the spectra of Lyapunov exponents, etc. 
We make comparison of the main concept of our paper with Turing's morphogenesis and John von Neumann 
automata, considering that this may be not only formal one, but will give ideas for the chaos development 
in the morphogenesis of Turing and for self-replicating machines. 

To provide rigorous study of the subject, we introduce new definitions such as chaotic sets of functions, 
the generator and replicator of chaos, and precise description of ingredients for Devaney and Li-Yorke chaos in 
continuous dynamics. Appropriate simulations which illustrate the morphogenesis phenomenon are provided. 

Keywords: Morphogenesis, Self-replication, Hyperbolic set of functions. Chaotic set of functions. Period- 
doubling cascade, Devaney's chaos, Li-Yorkc chaos, Intcrmittcncy, Chaotic attractor. Chaos control. The 
double-scroll Chua's attractor, Quasipcriodicity 



There are indicating properties of different types of chaos (sometimes they are called in- 
gredients). Morphogenesis, which is understood as replication of a prior chaos to the large 
community of systems, saving the indicating properties is proposed and mathematically ap- 
proved. 



1 Introduction 

The study of morphogenesis is one of the oldest of all the sciences, dating back to ancient Greece where 
Aristotle described the broad features of morphogenesis in birds, fish, and ccphalopods. Even at this early 
stage of scientific thinking, he understood that an animal's egg had the "potential" for its final form, but that 
it did not contain a miniature version of the form itself [2] . The structure of morphogenesis in an organism 
is presented in Figure [TJ 

In our paper, we try to follow the mechanism of morphogenesis, accepted as in classical and original 
meaning in biology [3]-[5l, not only in the sense of just replication of chaos of given type, but also saving the 
main assumption of Aristotle considering complexification of the form during its extension. That is, we try 
to establish the similarity between development from embryos to adult organisms considering development 
of chaos from simple observable systems to more complex ones. In the same time, we try to use the term 
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An embryo 



An adult organism 



Fi gur6 1 '. Morphogenesis of an organism 

morphogenesis not considering Aristotle's idea strictly, but issuing from the sense of words morph meaning 
"form" and genesis meaning "creation" [2]. In other words, similar to the ideas of Rene Thom [1], we employ 
the word morphogenesis as its etymology indicates, to denote processes creating forms. One of the purposes 
of our paper is to consider just the self-replication of chaotic systems without an emergence point of chaos, 
that is, we consider just the beginning of a technical task of formation of chaos of a given type through 
reproduction not necessarily considering complexification. This is different from Aristotle's idea in the sense 
that morphogenesis is effectively used in fields such as urban studies [B], architecture [7], mechanics [S], 
computer science [S] , linguistics [T^ and sociology [TTJ [T^ • 

The way of morphogenesis of chaos by the usage of a prime chaos as an embryo is presented in Figure [21 
We understand the prime chaos as that one which rigorously examined for the presence of chaos of a certain 
type. For instance a chaos of the logistic map obtained through period-doubling cascade can be considered 
as a prime chaos, and a system which provides the chaos will be called the generator (system). Similarly, 
the system which shares the type of chaos with the generator being influenced will be called the replicator 
(system). In the paper, we consider the replicator as non-chaotic and the sense of generator and replicator 
systems will be explained in the next section. 



A prime chaos 



A result of morphogenesis 



Figure 2: Morphogenesis of chaos 

Morphogenesis was deeply involved in mathematical discussions through Turing's investigations [M] as 
well as the relation of the concept of structural stability [13] . The paper [M] was one of the first papers that 
considers mathematically the self-replicating forms using a set of reaction-diffusion equations [15j . 

Morphogenesis of chaos as understood in our paper, is a form-generating mechanism emerging from a 
dynamical process which is based on self-replication of chaos. Here, we accept the form (morph) not only 
as a type of chaos, but also accompanying concepts as the structure of the chaotic attractor, its fractal 
dimension, form of the bifurcation diagram, the spectra of Lyapunov exponents, inheritance of intermittency, 
etc. Not all of the properties have been discussed in our paper, but wc suppose that they will be under 
investigation in our next papers. 

The concept of self-replicating machines, in the abstract sense, starts with the ideas of von Neumann [TB] 
and these ideas are supposed to be the origins of cellular automata theory [15] . According to von Neumann, 
it is feasible in principle to create a self-replicating machine, which he refers as an "automaton" , by starting 
with a machine A, which has the ability to construct any other machine once it is furnished with a set of 
instructions, and then attaching to A another component B that can make a copy of any instruction supplied 
to it. Together with a third component labeled C, it is possible to create a machine, denoted by i?, with 
components A, B and C such that C is responsible to initiate A to construct a machine as described by 
instructions, then make B to create a copy of the instructions, and supply the copy of the instructions to 
the entire apparatus. The component C is referred as "control mechanism". It is the resulting machine i?', 
obtained by furnishing the machine R by instructions //j, that is capable of replicating itself. Multiple usage 
of the set of instructions /jj is crucial in the mechanism of self- replication. First, the instructions must be 
fulfilled by the machine A, then they must be copied by _B, and finally the copy must be attached to machine 
R to form the system R' once again [151 E] . 

Our theory of morphogenesis of chaos relates the ideas of von Neumann about sclf-replicating machines 
in the following sense. Initially, we take into account a system of differential equations (the generator) 
which plays the role of machine A as in the ideas of von Neumann, and we use this system to influence in 
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a unidirectional way, another system (the rephcator) in the role of machine _B, in such a manner that the 
replicator mimics the same ingredients of chaos furnished to the generator. In the present paper, we use 
such ingredients in the form of period-doubling cascade, Devaney's and Li-Yorke chaos. In conclusion, the 
generator system with the replicator counterpart together, that is the result-system, admits ingredients of 
the generator. In other words, a known type of chaos is self-replicated. 

To understand the concept of our paper better, let us consider morphogenesis of fractal structures [TTIITS]. 
It is important to say that Mandelbrot fractal structures exhibit the appearance of fractal hierarchy looking 
in, that is, a part is similar to the whole. Examples for this are the Julia sets [TH [20] and the Sierpinski 
carpet [21] . In our morphogenesis both directions, in and out, are present. Indeed, the fractal structure 
of the prior chaos has hierarchy looking in, and the structure for the result-system is obtained considering 
hierarchy looking out. that is, when the whole is similar to the part. 

Two different mechanisms of morphogenesis are considered in our paper. In the first one, the generator is 
a first link of a chain where all other systems are consecutive links. This mechanism is illustrated in Figure 
[3] In another way, the generator connected with all other systems directly such that it is rather core than 
a beginning element and this method is described in Figure [4l In our paper we do not discuss theoretically 
constraints on the dimension of the result-system, but under certain conditions it seems the dimension is not 
restricted for both mechanisms. However, this is definitely true for the core mechanism even with infinite 
dimension. We call the two ways as the chain and the core mechanisms, respectively. We will discuss and 
simulate the chain mechanism in our paper, mainly, since the core mechanism can be discussed very similarly. 
One can invent other mechanisms of the morphogenesis, for example, by considering "composition" of the 
two mechanisms proposed presently. 
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Figure 3: Morphogenesis of chaos through consecutive repHcations 




Figure 4: Morphogenesis of chaos from a prior chaos as a core 

Successive applications of the morphogenesis procedure resemble the episodes of embryogenesis |3], in 
the biological sense. That is, in each "morphogenetic step", the regeneration of the same type of chaos in 
an arbitrary higher dimensional system is obtained. Moreover, elements inside the chaotic attractor of the 
generator, presented in the functional sense in our paper, play a role similar to the form-generating chemical 
substances, called morphogens, which are responsible for the morphogenetic processes |14| . In other words, 
these elements are responsible for the morphogenesis of chaos. 

In our paper, we are describing a process involving the replication of chaos which does not occur in the 
course of time, but instead an instantaneous one. In other words, the prior chaos is mimicked in all existing 
replicators such that the generating mechanism works through arranging connections between systems not 
with the lapse of time. 

Symbolic dynamics, whose earliest examples were constructed by Hadamard [22j and Morse [23], is one 
of the oldest techniques for the study of chaos. Symbolic dynamical systems are systems whose phase space 
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consists of one-sided or two-sided infinite sequences of symbols chosen from a finite alphabet. Such dynamics 
arises in a variety of situations such as in horseshoe maps and the logistic map. The set of allowed sequences 
is invariant under the shift map, which is the most important ingredient in symbolic dynamics |24j-|28|. 
Moreover, it is known that the symbolic dynamics admits the chaos in the sense of both Dcvancy and Li- 
Yorke [26l [29l |30l |31] |32] . The Smale Horseshoe map is first studied by Smale |33] and it is an example of a 
diffeomorphism which is structurally stable and possesses a chaotic invariant set [25l [26l [27] . The horseshoe 
arises whenever one has transverse homoclinic orbits, as in the case of the Duffing equation People 
used the symbolic dynamics to discover chaos, but we suppose that it can serve as an "embryo" for the 
morphogenesis of chaos. 

The phenomenon of the form recognition for chaotic processes began in pioneering papers already. It 
was Lorenz |35| who discovered that the dynamics of an infinite dimensional system being reduced to three 
dimensional equation can be next analyzed in its chaotic appearances by application of the simple unimodal 
one dimensional map. Smale [33] explained that the geometry of the horseshoe map is underneath of the Van 
der Pol equation's complex dynamics investigated by Cartwright and Littlewood [36] and later by Levinson 
|37j . Nowadays, the Smale horseshoes with its chaotic dynamics, is one of the basic instruments when one 
tries to recognize a chaos in a process. Guckenheimer and Williams gave a geometric description of the flow 
of Lorenz attractor to show the structural stability of codimension 2 [38]. In addition to this, it was found 
out that the topology of the Lorenz attractor is considerably more complicated than the topology of the 
horseshoe [53]. Moreover, Levi used a geometric approach for a simplified version of the Van der pol equation 
to show existence of horseshoes embedded within the Van der Pol map and how the horseshoes fit in the 
phase plane |39| . In other words, all the mentioned results say about chaos recognition, by reducing complex 
behavior to the structure with recognizable chaos. In [301 IMl HOI HB HH 1131 SI]; we provide a different and 
constructive way when a recognized chaos can be extended saving the form of chaos to a multidimensional 
system. In the present study we generalize the idea to the morphogenesis of chaos. 

Now, to illustrate the concept discussed in our paper, let us visualize the morphogenesis of chaos through 
simulating the Poincare sections of both generator and replicator systems. For our purposes we shall take 
into account the Dufhng's chaotic oscillator [45], represented by the differential equation 

x" + 0.05x' + = 7. 5cos{t). (1.1) 

Defining the variables xi = x and X2 = x' , equation ()1.1|) can be reduced to the system 

X'l = X2 2\ 

x'2 = -0.05x2 - A + 7.5cos(<). ^ • ' 

Next, we attach three replicator systems, by the chain mechanism shown in Figure [S] to the generator (|1.2p 
and take into account the 8— dimensional system 

x'l = X2 

x'^ ^ -0.05x2 - xf + l.hcos{t) 

X3 = X4 -I- Xi 

x\ = -3x3 - 2x4 - 0.008x| + X2 ,^ . 

x'^ = X6+X3 

x'e = -3x5 - 2.1x6 - 0.007x1 + X4 

x'j = Xs+ X5 

Xg = -3.1x7 - 2.2x8 - 0.006xf + xg. 

We suppose that system (|1.3p admits a chaotic attractor in the 8— dimensional phase space. By marking 
the trajectory of this system with the initial data Xi(0) = 2, X2(0) = 3, X3(0) = X5(0) = X7(0) = —1, X4(0) = 
Xq(0) = xg(0) = 1 stroboscopically at times that are integer multiples of 27r, we obtain the Poincare section 
inside the 8— dimensional space and in Figure [SJ where the morphogenesis is apparent, we illustrate the 
2— dimensional projections of the whole Poincare section. Figure[5ja) represents the projection of the Poincare 
section on the xi — X2 plane, and we note that this projection is in fact the strange attractor of the original 
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generator system p.2p . On the other hand, the projections on the planes x^ — Xi, x^ — xq, and x-z — xg presented 
in Figure [51^&), Figure [5Jc) and Figure [SJd), respectively, are the attractors of the replicator systems. 

One can see that the attractors indicated in Figure O (b) — (d) repeated the structure of the attractor 
shown in Figure EJa) , and the similarity between these pictures is a manifestation of the morphogenesis of 
chaos. 



a b 
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Figure 5: The picture in (a) not only represents the projection of the whole attractor on the xi — X2 plane, but also the 
strange attractor of the generator. In a similar way, the pictures introduced in (b), (c) and (d) represent the chaotic attractors 
of the first, second and the third replicators, respectively. The resemblance of the presented shapes of the attractors between 
the generator and the replicator systems reveal the morphogenesis of chaos. 

Next, let us illustrate in Figure [6l which also informs us about morphogenesis, the 3— dimensional pro- 
jections of the whole Poincare section on the X2 — X4 — and X3 — x^ — X7 spaces. One can see in Figure 
[SJa) and Figure [5]J6) the additional foldings which is not possible to observe in the classical strange attractor 
shown in Figure [5l (a). 

Despite we are restricted to make illustrations at most in 3— dimensional spaces, taking inspiration from 
FigureOand FigurelHl one can imagine that the structure of the original Poincare section in the 8— dimensional 
space will be similar through its fractal structure, but more beautiful and impressive than its projections. 
From this point of view, we are not surprised since these facts have been proved theoretically in the paper. 

To explain the morphogenesis procedure of our paper, let us give the following information. It is known 
that if one considers the evolution equation u' = L[u\ + Fit, u), where L[u\ is a linear operator with spectra 
placed in the left half of the complex plane, then a function Fit) being considered as input with a certain 
property (boundedness, periodicity, almost periodicity) produces through the equation the output, a solution 
with a similar property, boundedness/periodicity/almost periodicity. In particular, in our paper, we solved 
a similar problem when the linear system has eigenvalues with negative real parts and input is considered 
as a chaotic set of functions with a known type. Our results are different in the sense that the input and 
the output are not single functions, but a collection of functions. In other words, we prove that both the 
input and the output are chaos of the same type for the discussed equation. The way of our investigation 
is arranged in the well accepted traditional mathematical fashion, but with a new and a more complex way 
of arranging of the connections between the input and the output. The same is true for the control results 
discussed in the paper. 
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Fi gure 6l In (a) and (6) projections of the result chaotic attractor on X2 — X4 — xq and 3:3 — x^ — 3:7 spaces are respectively 
presented. One can see in (a) and (b) the additional foldings which are not possible to observe in the 2— dimensional picture 
of the prior classical chaos shown in Figure [S] (a). In the same time, the shape of the original attractor is seen in the resulting 
chaos. The illustrations in (a) and (b) repeat the structure of the attractor of the generator and the similarity between these 
pictures is a manifestation of the morphogenesis of chaos. 



One of the usage area of master-slave systems is the study of synchronization of chaotic systems [IHl UT] . 
In 1990, Pecora and Carroh [47| reahzed that two identical chaotic systems can be synchronized under 
appropriate unidirectional coupling schemes. Let us consider the system 

x' = G{x), (1.4) 

as the master, where x E M'', such that the steady evolution of the system occurs in a chaotic attractor. Now, 
let us take into consideration the slave system whose dynamics is governed by the equation 

v' = H{x,y), (1.5) 

when the unidirectional drive is established, with H{x,y) verifying the condition 

Hix,y) = Gix), (1.6) 

for y = X, and takes the form 

y' = G{y), (1.7) 

which is a copy of system (jl.4p . in the absence of driving. In such couplings, the signals of the master system 
acts on the slave system, but the converse is not true. Moreover, this action becomes null when the two systems 
follow identical trajectories [JB] . The continuous control scheme [351 HI] and the method of replacement of 
variables [SOI \5T\ can be used to obtain couplings in the form of the system (|1.4p + (|1.5p . Synchronization of 
a slave system to a master system, under the condition (|1.6p . is known as identical synchronization and it 
occurs when there are sets of initial data ^ C M'' and ^ C M'' for the master and slave systems, respectively, 
such that the equation 

lim ||.T(t)-y(t)|| =0, (1.8) 
holds, where {x{t),y{t)) is a solution of system (|1.4p + p. 51) with initial data (x(0),y(0)) e x 
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When the slave system is defined different than the master eounterpart in the absenee of driving, one 
can consider generalized synchronization. In this case, the dimensions of the master and slave systems are 
possibly different, say d and r respectively, and the definition of identical synchronization given by equation 
(|1.8p is replaced by 

lim U[x{t)] - y{t)\\ = 0, (1.9) 

where '!/'[a;(t)] is a functional relation which determines the phase space trajectory y(t^ of the slave system 
from the trajectory x(t) of the master. Here, it is crucial that the functional 'il)\x(t)\ has to be the same for 
any initial data (a;(0), y(0)) chosen from the set ^ x where the sets ^ C M'' and '¥ <ZW are defined in 
a similar way as in the case of identical synchronization. Generalized synchronization includes the identical 
synchronization as a particular case, and when is identity, one attains the equation (|1.8p . 

In our theoretical results of morphogenesis of chaos, we use coupled systems in which the generator system 
influences the replicator in a unidirectional way, that is the generator affects the behavior of the replicator, 
but not the converse. The possibility of making use of nonidentical systems in the replication of a known 
type of chaos furnished to the generator is an advantage of the procedure. On the other hand, contrary to 
the method that we present, in the synchronization of chaotic systems, one does not consider the type of the 
chaos that the master and slave systems admit. The problem that whether the synchronization of systems 
implies the same type of chaos for both master and slave has not been taken into account yet. 

In our results, we do not consider the synchronization problem as well the asymptotic properties, such 
as equations (|1.8p or (|1.9p . for the systems; but we say that the type of the chaos is kept invariant in the 
procedure. That is why the classes which can be considered with respect to this invariance is expectedly 
wider then those investigated for synchronization of chaos. Since there is no strong relation and accordance 
between the solutions of the generator and the replicator in the asymptotic point of view, the terms master 
and slave as well as drive and response are not preferred to be used for the analyzed systems. 

Nowadays, one can consider development of a multidimensional chaos from a low-dimensional one in 
different ways. One of them is chaotic itinerancy |52)-|59]. The itinerant motion among varieties of ordered 
states through high-dimensional chaotic motion can be observed and this behavior is named as chaotic 
itinerancy. In other words, chaotic itinerancy is a universal dynamics in high-dimensional dynamical systems, 
showing itinerant motion among varieties of low-dimensional ordered states through high-dimensional chaos. 
This phenomenon occurs in different real world processes: optical turbulence [54]; globally coupled chaotic 
systems [57l [58] ; non-equilibrium neural networks [55| [5^ ; analysis of brain activities [60] , ecological systems 
[M] . One can see that in its degenerated form chaotic itinerancy relates to intermittency [HH [53] , since they 
both represent dynamical interchange of irregularity and regularity. 

Likewise the itinerant chaos observed in brain activities, we have low-dimensional chaos in the subsystems 
considered and high-dimensional chaos is obtained when one considers all subsystems as a whole. The 
main difference compared to our technique is in the elapsed time for the occurrence of the process. In our 
discussions, no itinerant motion is observable and all resultant chaotic subsystems process simultaneously, 
whereas the low-dimensional chaotic motions take place as time elapses in the case of chaotic itinerancy. 
Knowledge of the type of chaos is another difference between chaotic itinerancy and our procedure. Possibly 
the present way of replication of chaos will give a light to the solutions of problems about extension of irregular 
behavior (crises, collapses, etc.) in interrelated or multiple connected systems which can arise in problems 
of classical mechanics [63] and electrical systems [Ml [65] , economic theory [66] , brain activity investigations 

In systems whose dimension is at least four, it is possible to observe chaotic attractors with at least two 
positive Lyapunov exponents and such systems are called hyperchaotic j67j . An example of a four dimensional 
hyperchaotic system is discovered by Rossler [68]. Combining two or more chaotic, not necessarily identical, 
systems is a way of achieving hyperchaos |69j-|71j. However, in the present paper, we take into account 
exactly one chaotic system with a known type of chaos, and use this system as the generator to reproduce the 
same type of chaos in other systems. On the other hand, the crucial phenomena in the hyperchaotic systems 
is the existence of two or more positive Lyapunov exponents and the type of chaos is not taken into account. 
In our way of morphogenesis, the critical situation is rather the replication of a known type of chaos. 
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2 Preliminaries 



Throughout the paper, the generator will be considered as a system of the form 

x' = F{t,x), (2.10) 

where F : M x M™ — >■ R™, is a continuous function in all its arguments, and the replicator is assumed to have 
the form 

y' = Ay + g{.,y), (2.11) 

where g : M'" x M" — > M", is a continuous function in all its arguments, and the constant n x n real valued 
matrix A has real parts of eigenvalues all negative. To describe better how the replicator, (|2.1ip . is involved 
in the morphogenesis process, we use the dot, symbol in the system. To insert the replicator in the 
morphogenesis process we replace the dot, by a; or a dependent variable of a replicator-predecessor and 
obtain the system of the form 

y' = Ay + g{x,y). (2.12) 

Let us remind that in the morphogenesis process of our paper one generator and several replicators may 
be involved, but in future one can consider process of morphogenesis where several generators are used. 

In general, in our paper, the replicator system (j2.1ip docs not admit a chaos. That is, we say that the 
replicator is non-chaotic. We understand the phrase "non-chaotic" in the sense that if "." is, for example, 
replaced by a periodic function, then system (j2.1ip does not produce a chaos, but admits a periodic solution 
and it is stable. 

We call the system which unites the generator system and several replicators, of type (|2.12p . in morpho- 
genesis process as the result-system of morphogenesis or in what follows just as the result- system. 

In what follows R and N denote the sets of real numbers and natural numbers, respectively, and the 

uniform norm ||r|| = sup ||ri'|| for matrices is used in the paper. 

||'i'||=i 

It is easy to verify the existence of positive real numbers N and uj such that ||e"^*|| < iVe~"*, t > 0, and 
these numbers will be used in the last assumption below. 

The following assumptions are needed throughout the paper: 

(Al) There exists a positive real number T such that the function F(t,x) satisfies the periodicity condition 

F{t + T,x) = F{t, x), for alH e R, a; e R™; 

(A2) There exists a positive real number Lq such that ||F(t,a:;i) — F{t,X2)\\ < ^ X2\\ , for all t G 

R,a;i,2:2 G M™; 

(A3) There exist a positive number Hq < oo such that sup |li^(t, j:)j| = Hq; 

(A4) There exists a positive real number Li such that \\g{xi,y) — g{x2,y)\\ > Li \\xi — X2\\ , for all xi,X2 G 
R™, y e K"; 

(A5) There exist positive real numbers L2 and L3 such that j|.g(a;i,?;) — g{x2,y)\\ < L2 \\xi — X2\\ , for all 
xi,X2 e R", y e M", and \\g{x,yi) - gix,y2)\\ < L3 \\yi - 2/2II , for aU x € R", yi,y2 G R"; 

(A6) There exist a positive number AIq < 00 such that sup ||g(x,?/)|| — Mq; 
(A7) NL3 - cj < 0. 

Remark 2.1 The results presented in the remaining parts are also true even if we replace the non- autonomous 
system (|2.10p by the autonomous equation 

x' = F{x), (2.13) 
where F : M™ — > M™, to > 3, is continuous, and with conditions which are counterparts of {A2) and {A3). 
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Using the theory of quasihnear equations |72| . one can verify that for a given solution x{t) of system 
(|2.10p . there exists a unique bounded on M solution y{t) of the system y' = Ay + g{x{t),y), denoted by 
y{t) = 4>x{t){t) throughout the paper, which satisfies the integral equation 



y{t) = / e 



:^^'-^^g{x{s),y{s))ds. 



(2.14) 



J —oo 



Our main assumption is the existence of a nonempty set of all solutions of system (|2.10p . uniformly 
bounded on K. That is, there exists a positive real number H such that sup ||a;(i)|| < H, for all x{t) G s/^- 



We note that for aU y{t) G s^y one has sup \\y{t)\\ < M, where M = 

Next, we reveal that if the set ^4 is an attractor with basin ^a,-, that is, for each x{t) G there exists 
x{t) G £/x such that \\x{t) — x{t)\\ as t oo, then the set si/y is also an attractor in the same sense. In 
the next lemma we specify the basin of attraction of £/y. Define the set 

^y — {y{t) I y{t) is a solution of the system y' = Ay + g{x{t), y) for some x{t) G } . (2-17) 



Lemma 2.1 'Wy is a basin of £/y. 

Proof. Fix an arbitrary positive real number e and let y{t) G '^y be a given solution of the system 
y' = Ay + g{x{t),y) for some x{t) G '^^x- In this case there exists x{t) G £/x such that \\x{t) — x{t)\\ — > 
as t -> oo. Let a = ^Jpfi^^^jij^^ and y{t) = (j>x(t){i)- Condition {AT) implies that the number a is posi- 
tive. Under the circumstances, one can find i?o(e) > such that \i t > Rq then \\x{t) — x{t)\\ < ae and 
N \\y{Ro) - y{Ro)\\ e(^^3-cj)t < q,^^ xhe functions y{t) and y{t) satisfy the relations 



Let us introduce the sets of functions 




(2.15) 



and 



= {ix{t),cl)xit){t)) I x{t)e£/x} 



(2.16) 




and 




respectively. Making use of these relations one has 



2/(t)-y(i) = e^(*-«°)(y(i?o)-y(i?o)) 




e^(*-««)(2/(i?o)-¥(i?o)) 
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Therefore, we have 



+ 



,A{t-Ro) 



URo)~y{Ro)\\ 



f 




\\9{x{s 




- g{x{s),y{s))\\ ds 












f 


pAit-s) 






g{x{s),y{s))\\ ds 



<Ne--<^'-^°^\\yiRo)-y{Ro)\\ 

+ I Ne~^^'-'^L:,\\y{s)^y{s)\\ds 

J Bo 

+ [ iVe-"(*-^)L2 \\x{s) - x{s)\\ ds 



+ 



NL2at 



-ut ( ut _ ujRo 



) 



-NL3 I e 
'Ro 



Let u : [i?o,oo) [0,oo) be a function defined as u{t) = e"* \\y{t) — y{t)\\ . Then we reach the inequahty 



(i)<A^e"^«||y(i?o)-y(i?o)|| 



NL2ae 



NL3 I u{s)ds. 
'Ro 



Now let = i^^^e"* and = ilit) + c where c = A^e"-"" \\y{Ro) -y{Ro)\\ - i^LLaSle"-"". Using 
these functions we get 



u{t) < (j){t) + NL3 / u{s)ds. 

JRo 

Applying Lemma (Gronwall) [73] to the last inequality for t > Rq, we obtain 

u{t) <c + iPit) + NL3 f e^^^^'^'^cds + NL3 I e^^-^(*-'V(s)(is 
Jro JRo 



and hence, 



u{t) < c + m + c (e^^^(*-^-) - 1) + f 1 _ piNLs-^Kt-Ro) ^ 



NL2ae 



NL2ae 
ui 



uj(lu — NL3) 

N\\y{Ro)-yiRo)\\e'^^°e^^^'-'-''"^ 

{NL3-iu){t-Ro')^ 



^uRo^NL^it-Ro) ^ L2L3ae I _^ _ ^ 

— NL3) 



Thus, 



hit) Vm < + N \\y{Ro) y(i?o)|| e(^^-"'(*-i^o) 



NL2ae 



,{NL3-u){t-Ro) _^ 



N'^L2L3ae 

u}(lu — NL3) 
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1 - e 



{NL3-co){t-Ro) 



NL2ae ^(jvL3-a;)(t-i?,o)^ N^L2L3ae _ ^{NL3-uj)(t-Ro) 



= N ||y(i?o) - y(i?o)|| e(^^3--)(t-i?o) + ^L^ae^ _ ^^^NL,-u)i,t-Ro)\^ 
< N ||y(i?o) - y(i?o)|| e(^^-^-)(*-^») + 



^3 

NL2ae 



Consequently, for < > 2Rq, we have — y{i)\\ < ^1 + ^j^^^^^ j ae = e, and hence \\y{i) — y{t)\\ as 
f — >■ 00. 

The proof is completed. □ 
Now, let us define the set 

^ = {{x{t), y{t)) I {x{t), y{t)) is a solution of system ((^1^ + ((TT^ such that a;(i) G '^^} . 

Next, we state the following corollary of Lemma [^TTI 

Corollary 2.1 '% is a basin of . 

Proof. Let {x{t), y{t)) G be a given solution of system (|2.10p + (|2.12|) . According to Lemma \2A\ one can 
find {x{t),y{t)) € such that \\x{t) — x{t)\\ — )■ as i — J- 00 and \\y{t) — y{t)\\ as t ^ 00. Consequently, 
\\{x{t),y{t)) - {x{t),y{t))\\ a.s t —?' 00. The proof is finalized. □ 



3 Description of chaotic sets of functions 

In this section of the paper, the descriptions for the chaotic sets of continuous functions will be introduced 
and the definitions of the chaotic features will be presented, both in the Devaney's sense and in the Li-Yorke 
sense. 

Let us denote by 

^ = {'ijj{t) I -0 : R ^ /\ is continuous} (3.18) 

a collection of functions, where if C M^, g G N, is a bounded region. 

We start with the description of chaotic set of functions in Devaney's sense and next continue with the 
Li-Yorke counterpart. 



3.1 Chaotic set of functions in Devaney's sense 

In this part, we shall elucidate the ingredients of the chaos in Devaney's sense for the set ^ and the first 
definition is about the sensitivity of chaotic set of functions. 

Definition 3.1 £^ is called sensitive if there exist positive numbers e and A such that for every i/jlt) £ ^ 
and for arbitrary S > there exist ipit) € to (zM. and an interval J d [to, 00) with length not less than A 
such that \\ip{to) — "0(^0)11 < S and ~ V'(0|| ^ '^^^ t £ J. 

Definition 13.11 considers the inequality (> e) over the interval J. In the Devaney's chaos definition for the 
map, the inequality is assumed for discrete moments. Let us reveal how one can extend the inequality from 
a discrete point to an interval by considering continuous flows. 

In [26], it is indicated that a continuous map (p : A A, with an invariant domain A C M'^ , A: e N, has 
sensitive dependence on initial conditions if there exists e > such that for any x € A and any neighborhood 
^ of X, there exist y G ^ and a natural number n such that ||(/3"(a;) — Lp"{y)\\ > e. 

Suppose that the set x/x satisfies the definition of sensitivity in the following sense. There exists e > 
such that for every x(t) g jz/x and arbitrary S > 0, there exist x{t) G £/x, to G R and a real number C ^ 
such that ||x(to) — ^(^o)|| < S and ||a;(C) — a^(C)ll > £■ 
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In this case, for given x{t) £ J24 and 5 > 0, one can find x{t) G s^^. and C > such that |la;(to) ^ ^(^o)ll < ^ 
and — x(C)|| > e. Let A = g^j— and take a number Ai such that A < Ai < jjf^;- Using the appropriate 

integral equations, for t € [C, C + we have 



[F{s,x{s))~F{s,x{s))]ds 



C 



\\x{t)-xm>MO-m\\- 

>e-2HLoAi 
e 

> -. 

- 2 

The last inequality confirms that satisfies the Definition 13.11 with e = | and J = [Cj C + ^i]- So the 
definition is a natural one. It provides more information then discrete moments and for us it is important 
that the extension on the interval is useful to prove the property for morphogenesis process. 

In the next two definitions, we continue with the existence of a dense function in the set of chaotic 
functions followed by the transitivity property. 

Definition 3.2 £i§ possesses a dense function i}}* {t) G if for every function 'il'{t) € ^M, arbitrary small 
e > and arbitrary large E > 0, there exist a number ^ > and an interval J C M ujith length E such that 
\\4it) -ij*{t + C)\\ < e, for allte J. 

Definition 3.3 £i§ is called transitive if it possesses a dense function. 

Now, let us recall the definition of transitivity for maps |26j . A continuous map (p with an invariant 
domain A C M'^, fc G N, possesses a dense orbit if there exists c* € A such that for each c G A, arbitrary small 
real number e > and arbitrary large natural number Q > 0, there exist natural numbers kg and such 
that ||(y9"(c) — (^"+'^°(c*)|| < e, for each integer n between and + Q, and maps which have dense orbits 
are called transitive. 

Suppose that satisfies the transitivity property in the following sense. There exists x* (t) G ^ such 
that for each x{t) G ^4, arbitrary small positive real number e and arbitrary large natural number Q, there 
exist natural numbers fco and Iq such that \\x{nT) — x*{{n + fco)r)|| < e, for each integer n between and 
lo + Q. 

Fix arbitrary x{t) G £/x, e > and Q G N. Under the circumstances, one can find fco,^o G N such that 
||a;(nT) - x*{{n + ko)T)\\ < ee'^^^'^, for each integer n = Iq, . . . Jo + Q- 

Using the condition {A2) together with the convenient integral equations that x{t) and x*{t) satisfy, it is 
easy to obtain for t G [IqT, {Iq + Q)T] that 

\\x{t) - x*{t + koT)\\ < \\x{loT) - x*{{lo + ko)T)\\ + f Lo \\x{s) -x*{s + koT)\\ ds, 

JloT 

and by the help of the Gronwall-BcUman inequality [74] , we get 

\\x{t) -x*{t + koT)\\ < \\x{loT) - x*{{lo + ko)T)\\ e^o(*-'o^) < e. 



The last inequality shows that the set ^/^ satisfies the Definition 13.21 with ^ = koT and E = QT, and is 
transitive according to Definition 13.31 



Definition 3.4 admits a dense countable collection G of periodic functions if for every function 
ipit) G ^1 arbitrary small e > and arbitrary large _B > 0, there exist ipit) G o,nd an interval J C M, with 



length E, such that 



m - m 



< e, for all t G J. 



Let us remind the definition of density of periodic orbits for maps [26] . The set of periodic orbits of a 
continuous map (p with an invariant domain A C M*^, fc G N, is called dense in A if for each c G A, arbitrary 
small positive real number e and arbitrary large natural number Q there exist a natural number Iq and a 
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point c G A such that the sequence {(y9*(c)} is periodic and |lv"(c) — <y5"(c)|| < e, for each integer n between 
Iq and Iq + Q. 

Let us denote by the set of all periodic functions inside £/x- Suppose that satisfies density of 
periodic solutions as follows. For each function x{t) G £/x, arbitrary small e > and arbitrary large Q G N, 
there exist x{t) e and a number Iq (^N such that \\x{nT) — ■r(r7,r)j| < e, for each integer n between Iq and 
lo + Q- 

Let an arbitrary function x{t) S £/x, arbitrary small e > and arbitrary large Q G N be given. In that 
case, there exist fco, /q € N such that ||a;(r7,r) — ■r(nT)|| < ee~^°'^^ , for each integer n = Iq, ■ ■ ■ Jo + Q- 
It can be easily verified that for t £ [loT, {Iq + Q)T] the inequality 

\\xit) ~ x{t)\\ < \\x{loT) - x{loT)\\ + f Lo \\xis) - x{s)\\ ds, 

JloT 

holds and therefore we have 

\\x{t) - x{t)\\ < \\x{loT) - x{loT)\\ eio(*-'o^) < e. 
Consequently, the set ^/^ satisfies Definition 13.41 with E = QT and J = [IqT, IqT + E]. 
Definition 3.5 The collection of functions SS is called a Devaney's chaotic set if 
(Dl) ^ is sensitive; 
(D2) ^ is transitive; 

(D3) ^ admits a dense countable collection of periodic functions. 

In the next discussion, the chaotic properties on the set M will be imposed in Li-Yorke sense. 

3.2 Chaotic set of functions in Li-Yorke sense 

The ingredients of Li-Yorke chaos for the collection of functions which is defined by (|3.18p . will be described 
in this part of the paper. Making use of discussions similar to the ones made in the previous subsection, we 
extend, below, the definitions for the ingredients of Li-Yorke chaos from maps [ZS]-IZH] to continuous fiows 
and we just omit these indications here. 

Definition 3.6 A couple of functions ■0(0) & ^ 'X ^ is called proximal if for arbitrary e > 0, E > 

there exist infinitely many disjoint intervals of length not less than E such that — V'lOU < /'^'^ each t 

from these intervals. 

Definition 3.7 A couple of functions {il>{t),il>{t)) € ^ is frequently [e^, IS.) ^ separated if there exist pos- 
itive numbers eg, A and infinitely many disjoint intervals of length not less than A, such that — ''P{t)\\ > 
eq /o'' each t from these intervals. 

Remark 3.1 The numbers Eq o,'>T'd A taken into account in Definition \3.7\ depend on the functions ip(t) and 

m- 

Definition 3.8 A couple of functions (ip{t),ilj{t)) £ ^ x ^ is a Li~Yorke pair if they are proximal and 
frequently {eo, A)— separated for some positive numbers eq o.iT'd A. 

Definition 3.9 An uncountable set <Z is called a scrambled set if '£ does not contain any periodic 
functions and each couple of different functions inside x 'lo is a Li— Yorke pair. 

Definition 3.10 ^ is called a Li— Yorke chaotic set if 

(LYl) There exists a positive real number Tq such that SS admits a periodic function of period kTo, for any 
fc e N; 

(LY2) £^ possesses a scrambled set 

(LY3) For any function ^{t) S ^ and any periodic function tl'{t) S l3§, the couple (^{t),^p{t)) is frequently 
(eo. A)— separated for some positive real numbers sq and A. 
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4 Hyperbolic set of functions 

The definitions of stable and unstable sets of hyperbolic periodic orbits of autonomous systems are given 
in [73] , and information about such sets of solutions of perturbed non-autonomous systems can be found in 
[50] , Moreover, homoclinic structures in almost periodic systems were studied in [51]- [55]. In this section, 
we give definition for an hyperbolic collection of uniformly bounded functions and before this, we start with 
definitions of stable and unstable sets of a function. 

We define the stable set of a function ^{t) £ SS, where the collection is defined by p.lSp . as the set of 
functions 



Definition 4.1 The set of functions SS is called hyperbolic if each function inside this set has nonempty 
stable and unstable sets. 

Theorem 4.1 If is hyperbolic, then the same is true for .s^/y. 

Proof. Fix arbitrary e > 0, y{t) = (j)x{t)it) € and let a = ^J^~l^+nl2 ' ^ = TTFrf - condition {A7), 
one can verify that the numbers a and /3 are both positive. 

Due to hyperbolicity of the function x{t) has a nonempty stable set W{x{t)) and a nonempty unstable 
set Wixit)). 

Take a function u{t) G W''{x{t)). Since \\x{t) — u{t)\\ — > as i — )• oo and NL^ — < 0, there exists a 
positive number which depends on e, such that \\x{t) — u{t)\\ < ae and e^^^^""'* < for t> Ri. 

Let us denote y{t) = (pu{t){i) and we shall prove that the function y{t) belongs to the stable set of y{t). 
The bounded on M functions y{t) and y{t) satisfy the relations 



W {^{t)) {u{t) e ^ I \\u{t) - i:{t)\\ ^ as t ^ cx)} , 



(4.19) 



and, similarly, we define the unstable set of a function tp{t) G ^ as the set of functions 



W {ij{t)) = {v{t) e ^ I \\v{t) - ^(t)\\ ^ as t ^ -oo} . 



(4.20) 




and 




respectively, for t>Ri. 

Therefore one can easily reach up the equation 



y{t)^y{t) = / e^<^'-%{x{s),y{s))-g{u{s),y{s))]ds 



J —OO 



+ f e^^*-^' {[g{x{s),y{s))-g{x{s),y{s))] + [g{x{s),y{s)) - g{u{s)M)]) ds. 



which implies that 
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Using the Gronwall type inequality indicated in |84| . we attain the inequahty 



For that reason, for t > 2i?i, one has 

< IH ae 

\ uj-NlJ 
= e. 

The last inequality reveals that \\y{t) — y{t)\\ as i — >■ 00, and hence y{t) belongs to the stable set 
W^yit)) ofy(i). 

On the other hand, let v{t) be a function inside the unstable set W'^{x{t)). 

Since \\x{t) — v{t)\\ — ?> as i — >■ —00, there exists a negative real number i?2(e) such that \\x{t) — v{t)\\ < 
/3e, for t < i?2- Let y{t) = (j)y(^t){t)- Now, our purpose is to show that y{t) belongs to the unstable set of y{t). 
By the help of the integral equations 



y{t)= e^<^'~^^gix{s),y{s))ds, 



and 



we obtain that 



—00 



t 

e' 

—00 
t 



yit) - m = / e^(*-^) [g{x{s), y{s)) - y{s))]ds 

J —00 

3^^'-^\g{x{s),y{s)) - g{v{s),y{s))]ds 

e^<^*-^\g(vis),y{s)) ~ g(v{s),yis))]ds. 
Therefore, for t < R2, one has 

\\y{t)-ym< f NL2e-'^^'-^^\\x{t)-v{t)\\ds 



t 



Hence, 



and 



< + sup \\y{t) - y{t)\\ 

W UJ t<R2 



„ , > ~. , n NLnBe NL-i ,, , , 
sup \\y{t) - y{t)\\ < — + sup \\y{t) - y{t)\\ 



sup \\y{t)-y{t)\\ < ^^frl <e- 
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The last inequality confirms that \\y{t) — y{t)\\ as f ^ — oo. Therefore y{t) € W'^{y{t)). 

Consequently, ^/y is hyperbolic since y{t) possesses both nonempty stable and unstable sets, denoted by 
W'^(y{t)) and W"{y{t)), respectively. 

The theorem is proved. □ 

Theorem 14.11 implies the following corollary. 

Corollary 4.1 If s^x is hyperbolic, then the same is true for . 

Next, we continue with another corollary of Theorem (|4.ip . following the definitions of homoclinic and 
heteroclinic functions. 

A function i!i{t) G ^ is said to be homoclinic to the function "00 (i) € i'o{t) 7^ if S 

On the other hand, a function V'2(i) € =^ is called heteroclinic to the functions ipo{t),il}i{t) G 3§, V'o(i) 7^ 

Corollary 4.2 If xi{t) € s^^ is homoclinic to the function Xa{t) € j^/x, XQ{t) 7^ x{t), then (pxi(t){i) G ■^y is 
homoclinic to the function 4'xo{t){i) € -c/j,. 

Similarly, if X2(t) € J^/^ is heteroclinic to the functions xo{t),xi{t) € £/x, Xo(t) 7^ X2{t), xi{t) 7^ 2^2(^)7 
then 4>x2{t)it) is heteroclinic to the functions (j)xo{t)it)i 4'xi(t){i) S ■^y 

5 Self-replication of chaos 

In this section, which is the main part of our paper, we theoretically prove that the set s^y replicates the 
ingredients of chaos furnished to the set J24 and as a consequence the same is also valid for the set . We 
start our discussions with chaos through period-doubling cascade, and continue with Devaney's and Li-Yorke's 
chaotic ingredients, respectively. 

The mechanism of morphogenesis was present before through the usage of the logistic map as an "embryo" 
to generate the chaotic structure embedded in a system with an arbitrary finite dimension. Devaney's chaos 
[30] , Li-Yorke chaos [43] , chaos through period-doubling cascade [44] and intermittcncy [30] were provided in 
a quasilinear system of the type 

Z'{t) ^ Az{t) + f{t, Z) + V{t, to, m) (r ni ^ 

z{to)^zo,{to,zo)e[0,l]xR^, ' 

where 



.it,to,f.) = ( < ' S^+iJ^'^l (5.22) 

^'"''^^ \ TOi, if C2,-i(io,Ai) < ^ < C2,:(io,/^), ^ ^ 

i is an integer, mo, mi € R", the sequence ({to,fi) = {Cj(^o,m)} is defined through equations Cj(^o,M) = 
i + Ki{tQ, fi), with Ki+i(to, n) = h{Ki{to, ji), noito, /i) = to, and h{x, /i) = fix{l — x) is the logistic map. In 
the present paper, we do not consider an "embryo" to generate the chaotic structure embedded in a system, 
but we make use of a chaotic system to replicate the same type of chaos. 

5.1 Period-doubling cascade and morphogenesis 

We start this subsection by describing the chaos through period-doubling cascade for the set of functions £/x 
and deal with its replication by the set of functions ^/y, which is defined by equation (|2.15p . 

Suppose that there exists a function G : M. x R™ x M — > which is continuous in all of its arguments 
such that G{t + T,x,fj,) = G{t,x,fj.), for alH g M, a; S W^,fj. e R, and F{t,x) = G{t,x,fj,oo) for some finite 
value /Ltoo of the parameter /x, which will be explained below. 
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To discuss chaos through pcriod-doubhng cascade, let us consider the system 



x' = G{t,x,n). (5.23) 

We say that the set s!/^ is chaotic through period-doubhng cascade if there exists a natural number k and 
a sequence of period-doubling bifurcation values {/im} , fJ-m — > Moo as m — > oo, such that for each m e N as 
the parameter fi increases (or decreases) through fim, system (|5.23|) undergoes a period-doubling bifurcation 
and a periodic solution with period fc2"T appears. As a consequence, at /i = ^oo, there exist infinitely many 
periodic solutions of system (|5.23p . and hence of system (|2.10p . all lying in a bounded region. In this case, 
the set £/x admits periodic functions of periods kT, 2kT, AkT, 8kT, ■ ■ ■ . 

Now, making use of the equation (j2.14p . one can show that for any natural number p, if x{t) e ^ is 
a pT— periodic function then (t>x{t)(t) € J^y is also pT— periodic. Moreover, condition {A4) implies that the 
converse is also true. Consequently, if the set admits periodic functions of periods kT, 2kT, AkT, 8kT, • • • , 
then the same is valid for £/y, with no additional periodic functions of any other period, and this provides us 
an opportunity to state the following theorem. 

Theorem 5.1 // the set ^ is chaotic through period- doubling cascade, then the same is true for siy. 

The following corollary states that the result-system (|2.10p + (|2.12p is chaotic through the period-doubling 
cascade, provided the system (|2.10p is. 

Corollary 5.1 // the set s^/x is chaotic through period- doubling cascade, then the same is true for si . 

Our theoretical results show that the existing replicator systems, likewise the generator counterpart, 
undergoes period-doubling bifurcations as the parameter /i increases or decreases through the values /im, m S 
N. That is, the sequence {/im} of bifurcation parameters is exactly the same for both generator and replicator 
systems. Since the generator system necessarily obeys the Feigenbaum universality i45i i85l [86l I87j . one can 
conclude that the same is also true for the replicator. In other words, when limm-j-oo ^"^^'"^^ is evaluated, 
the universal constant known as the Feigenbaum number 4.6692016 ... is achieved and this universal number 
is the same for both generator and replicators, and hence for the result-system. One can percieve this result 
as another aspect of morphogenesis of chaos. 

Next, let us illustrate by simulations, the morphogenesis of chaos through period-doubling cascade. In 
paper [88], it is indicated that the Buffing's equation 

x" + 0.3a;' + a;3 = 40cosi (5.24) 

admits the chaos through period-doubling cascade. Defining the new variables x\ = x and X2 = x' , equation 
(|5.24p can be rewritten as a system in the following form 



(5.25) 



2-1 = X2 

x'^ ^ -0.3x2 - arf -I- 40 cost. 
Making use of system (|5.25p as the generator, one can constitute the 8— dimensional result-system 

x'l ~ X2 

x'2 = -0.3a;2 ~ x\^ 40 cost 
4 ^ 2x3 - 2:4 + 0.4 tan (^^4^) 

x\ = I7.T3 - 6x4 -I- a;2 ,^ 

ccg = — 2x5 -I- 0.5sina:;6 — 4x4 ' 

Xg = -X5 - 4x6 - tan (^) 

x'7 = 2x7 + 5x8 - 0.0003(x7 - xg)^ - 1.5x6 

Xg = —5x7 — 8xg -t- 4x5, 

which involves three consecutive 2— dimensional replicator systems in addition to the generator, as the chain 
procedure indicated in Figure [3 
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Wc note that system (|5.26p exhibits a symmetry under the transformation 



J(f : (a;i,X2,a;3,X4,a;5,a;6,X7,a;8,i) (-Xi, -a;2, -Xs, -X4, -X5, -xg, -xy, -xg, t + tt). (5.27) 

According to our theoretical discussions, the result-system (|5.26p admits a chaotic attractor in the 
8— dimensional phase space which is obtained through period-doubling cascade. In Figure [3 (a) — (d), we 
illustrate the 2— dimensional projections of the trajectory of system (j5.26p . with the initial data a;i(0) = 
2.16, X2(0) = -9.28,0:3(0) = -0.21,0:4(0) = -2.03,2:5(0) = 3.36, X6(0) = -0.52,2:7(0) = 3.07,2:3(0) = -0.32, 
on the planes 2:1 — 2:2, 2:3 — 2:4, 2:5 — 2:6, and 2:7 — 2:8, respectively. The picture in Figure [71 (a), shows in fact 
the attractor of the prior chaos given by the system (|5.25p and similarly the illustrations in Figure [3 (6) — (d) 
correspond to the chaotic attractors of the first, second and the third replicator systems, respectively. The re- 
semblance between the shapes of the attractors of generator and replicator systems reflect the morphogenesis 
of chaos in the result-system (|5.26p . 



a b 




-6-4-2 2 4 6 -1.5 -1 -0.5 0.5 1 1.5 



d 




-5 5 -8 -6 -4 -2 2 4 6 8 



Fi gure 7: 2— dimensional projections of the chaotic attractor of the result-system I I5.26I I. (a) Projection on the x\ — X2 plane, 
(b) Projection on the X3 — X4 plane, (c) Projection on the X5 — xq plane, (d) Projection on the — xg plane. The picture in 
(a) shows the attractor of the prior chaos produced by the generator system 15.2511 and in (b) — (d) the chaotic attractors of the 
replicator systems are observable. The illustrations in (6) — (d) repeated the structure of the attractor shown in (a), and the 
similarity between these pictures is an indication of the morphogenesis of chaos. 

To obtain better impression about the chaotic attractor of system (|5.26p . in Figure |8] we demonstrate the 
3— dimensional projections of the trajectory with the same initial data as above, on the 2:3 — 2:5 — 2:7 and 
2:4 — 2:6 —2:8 spaces. Although we are restricted to make illustrations at most in 3— dimensional spaces and not 
able to provide a picture of the whole chaotic attractor, the results shown both in Figure [7] and Figure [8] give 
us an idea about the spectacular chaotic attractor of system ()5.26p . We note that the presented attractors 
are symmetric around the origin due to the symmetry of the result-system (|5.26p under the transformation 
^ introduced by ([07]) . 

Now, let us show that the first replicator system which is included inside (|5.26p satisfies the condition 

Let us denote by ||.|| the matrix norm which is induced by the usual Euclidean norm in M.K That is, 

||r|| = maxjy^ : <; is an eigenvalue of F^^F} (5.28) 
for any I x / matrix F with real entries, and F^ denotes the transpose of the matrix F |89j . 
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b 





Figure 8: 3— dimensional projections of the chaotic attractor of the result-system 115.261 1. (a) Projection on the X3 — — x-j 
space, (b) Projection on the X4 — xe — space. The illustrations presented in (a) and (b) give information about the impressive 
chaotic attractor in the 8— dimensional space. 



When the system 



4 = 2x3 - a;4 + 0.4 tan 



YJx-i — 6x4 



2^2 



(5.29) 



is considered in the form of system (j2.12p . one can sec that the matrix A can be written in the form 
2 -1 . 

which admits the complex conjugate eigenvalues —2 =p i. 



A 



17 -6 

The real Jordan form of the matrix A is given by J 
1 



satisfied where P = 



1 4 



-2 -1 
1 -2 

Evaluating the exponential matrix e'^* we obtain that 



and the identity P ^AP = J is 



cos{t) ~sin{t) 
sin(t) cos{t) 



P- 



(5.30) 



Taking N = \\P\\ ||P ^|| < 18 and w = 2, one can see that the inequality < Ne holds for all 

t > 0. The function 5 : x ^ defined as 



satisfies the conditions (AA) and {A5) with constants Li 



0.4 tan ( I ,X2 



10 



50 



, L2 = \/2 and L3 = 0.08 since the chaotic 



attractor of system (|5.26p satisfies the inequalities \xi \ < 6, \x3\ < |, and consequently |^-^^^| < f . Therefore, 
the condition {A7) is satisfied. 

In a similar way, for the second replicator system, making use of l^sj < | once again, one can show that 
the function /i : x ^ defined as 

h{x3, 2:4, X5, ccg) = 1^0.5 sinxe — 4a;4, — tan ) 
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satisfies the counterparts of the conditions (Ai) and {A5) with constants Li ~ ^2 = 4\/2 and L3 = i. 
On the other hand, when one considers the third replicator, it is easy to show that the function k : 
X M2 ^ ]r2 defined as 

/c(x5, xq,X7, xs) = (0.0003(2:7 — xg)^ — 1.5x6, 4x5) 

satisfies the counterparts of the conditions (A4) and {A5) with constants Li = L2 = 4-\/2 and L3 = 0.19 
since the chaotic attractor of system (|5.26p satisfies the inequalities \x7\ < 8, \xs\ < 4. 



5.2 Machinery of morphogenesis of Devaney's chaos 

In this part of our paper, we will prove theoretically that the ingredients of Devaney's chaos furnished to the 
set £/x are all replicated by the set £/y. 

Lemma 5.1 Sensitivity of the set six implies the same feature for the set s/y. 

Proof. Fix an arbitrary S > and let y{t) e s/y be a given solution of system ()2.12|) . In this case, there 
exists x{t) €z such that y(t) = (t>x(t){t)- 

Let us choose an e = l{d) > small enough which satisfies the inequality ^1 + Jl^\,^ ^ e < 5. Then take 
R = R{e) < sufficiently large in absolute value such that MML(.{'^-nL:,)r ^ ^^^^ ^et 5i = 5i{e, R) = ee^«^. 
Since the set of functions ^ is sensitive, there exist positive real numbers eo and A such that the inequalities 
||x(to) ~2;(to)|| < 5i and \\x{t) ~x{t)\\ > eo, t E J, hold for some x{t) £ ^4, io G and for some interval 
J C [to, 00) whose length is not less than A. 

Using the couple of integral equations 

x{t)=x{to)+f F{s,x{s))ds, 

J to 



x{t) = x{to)+ / F{s,x{s))ds 

J to 

together with condition {A2), one see that the inequality 



\\x{t)^x{t)\\ < \\x{to)-xito)\\ + 



Lq ||x(s) — x{s)\\ ds 



holds for t E [to + R, to] . Applying Gronwall-Bellman inequality [74] , we obtain 

\\xit)-x{t)\\ < \\x{to)-xito)\\e'^o\t-to\ 

and therefore \\x{t) — x{t)\\ < e, for t E [to + R,to]- 

Suppose that y{t) = (j>x(t){i) & First, we wiU show that ||2/(io) — y(to)|| < ^■ 
The functions y(i) and yit") satisfy the integral equations 



and 



g{x{s),y{s))ds 



g{x{s),y{s))ds, 



respectively. Therefore, 
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y{t) - v{t) - / e-4(*-^) [g{x{s), y{s)) - y(s))]ds 



and hence we obtain 



\\y{t)-ym< / Ne-'^<-'-^'^\\g{x{s),y{s))~g{xis),y{smds 



to+R 



+ / Ne--<^'-^^\\g{x{s),y{s))~g{x{s),y{s))\\ds 
Jto+R 
pto+R 

+ / ^e-(*-^) \\g{x{s),y{s)) - g{x{s),y{sm ds. 



Since \\x{t) — x{t)\\ < e for i e [to + ^, to], one has 



Mt) - < NL^ / e-'^(*-^) |ly(,s) - y{s)\\ ds 



to+R 

+NL2e [ e~"(*-")ds + 2MoN 

Jto+R 



to+R 



< NL3 
NL2e 



e ~- lly(s) - y{s)\\ ds 



to+R 



-Ljt^^ut _ ^u{to+R)\^ _^ '^^'^0^ ^-u{t-to-R)^ 



Now, introduce the functions u{t) = e"' -y{t)\\ , fc(t) = ii^e"*, and = c + where c 

U UJ J 

These definitions give us the inequahty 



u{t) < h{t) + / NL3u{s)ds 



to+R 



Applying Lemma 2.2 [90] to the last inequality, we get u{t) < h{t) + NL^} I 

Jt( 

Therefore, for i e [io + R, to\ we have 



^NL3{t-s) 



to + fl 



h{s)ds. 



Lit) < c + kit) + c (e^i3(t-to-fl) _ l^ + N^L^L-se ^^L^t I ~ ^(.-nl,)s^^ 

V J LO ./t„_LR 



kit) + ce^^^(*-*«"^' + J^!M^e"* 
w(ti; — NL3) 

NL2e /2MoN TVLse 



1 - e 



to+-R 
{NL3~u){t^to~R) 



V 



^uR^NLa{t~to~R) 



N'L2L3e 
f 

a;(a; — NL3) 



1 - e 



(NL3-u:)(t~to-R) 



Hence, 



\\yit)-yit)\\ < 

^ N^L2L3e 



NL2I f2MoN TVLae 



^iNL3-uj){t-to-R) 



1 - e 



{NL3-Lj){t-to-R.) 
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NL 



2e 



UJ - NL3. 

NL^e 



2MoN 



NL 
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CJ - NL3 

The last inequality leads to 



U! UJ — NL3 



^{NL3-ui){t^to~R) 



\\yit)^ym< 



UJ - NL3 



2MoN 



,{NL3-iu){t-to-B) 



Consequently, 



\y{to) - y(to)ll < ^ + — ^e("-^^3)fl 

UJ — I\ L3 U! 

NL2 



< 5. 



UJ - NL3 



In the remaining part of the proof, we will show the existence of a positive number ei and an interval 
C J, with a fixed length which is independent oi y{t),y{t) e £/y, such that the inequality \\y{t) — y{t)\\ > t\ 
holds for ah t e 

/ .91 (a;, y) \ 



Suppose that g(x^ y) 



.92 (a;, y) 



where each gj, 1 < j < n, is a real valued function. 



\ 9n{x,y) J 

Since for each x{t) E the functions x'{t) and '/'^(j) (0 inside the tubes with radii Hq and A/ + A/q, 
respectively, one can conclude that the set of functions jz^. and £/y are both equicontinuous on M. Making use 
of the uniform continuity of the function p : R"'xR™xR" M" defined as X2, X3) = g{xi,X3) — g{x2,X3) 
on the compact region 

^ = {(xi,a;2,a;3)eM"xM"xR" | ||a;i|| < H,\\x2\\ < H,\\x3\\ < M} , 

together with the equicontinuity of the sets ^/^ and x/y, one can easily verify that the set 

^ = {.g,(a;(t),0,(t)(i)) - W) I 1 < J < x{t) € € (5.31) 

is an equicontinuous family on R. 

Therefore, there exists a positive real number t < A, independent of the functions x{t),x[t) G £/x, 
y{t),y{t) G si/y, such that for any ^1,^2 G R with \ti — t2| < t the inequality 

Lieo 



1(5, (x(ti),y(ti)) -5, (x(ti),2/(ti))) - (.gj (.T(t2), 2/(^2)) - 5, (^(ts), 2/(^2)))! < 



2n 



(5.32) 



holds, for all 1 < J < n. 

Condition (AA) implies that, for all f G J, the inequality \\g(x{t),y{t)) — g{x{t),y{t))\\ > Li \\x{t) — x{t)\\ 
holds. Therefore, for each t G J, there exists an integer jo = joit), 1 < jo < n, such that 

\9Mt).yit))-9jom,ym > ^ \\xit)-xit)\\ . 

Otherwise, if there exists s G J such that for all 1 < j < n, the inequality 

\gjix{s),y{s))~g,{x{s),y{s))\ < ^ ||x(s) - a;(s)|| 
holds, then one encounters with a contradiction since 

\\g{x{s),y{s)) - g(x{s),y{s))\\ 

n 

< \9ji^is),yis)) - gj{x{s),y{s))\ 

< Li \\x{s) -x{s)\\ . 
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Now, let So be the midpoint of the interval J and 9 — sq — ^. One can find an integer jo = jo('So)i 
1 !i jo < such that 

l5jo(a;(so),2/(so))-5,o(^(5o),y(so))|> ^||x(so)-x(so)|| > (5.33) 
On the other hand, making use of inequality (|5.32p . for all t G [6*, + t] we have 

ISjo (a;(so),2/(so)) - 5io (^(so), y(so))| - Isjo {■•>:{t),y{t)) ~ 9jo (^W, 2/(0)1 

< Kffjo {x{t)^y{t)) - 930 ix(t),y{t))) - (ffio (2;(so),2/(so)) -5io {x{so),y{so)))\ 



< 



2n 



and therefore by means of (|5.33[) . we obtain that the inequality 

\9jo ixit),y{t)) - gjo {xit),y{t))\ > \gj^ {x{so),y{so)) - g^^ (x(so), 2/(so))| - 



2n ~ 2n 



(5.34) 



holds for ah t e [6*, 6* + r] . 

By applying the mean value theorem for integrals, one can find Si, S2, . . . , s„ G [9,0 + t] such that 



[9{x{s),y{s)) - g{x{s), y(s))] ds 



[gi{x{s),y{s)) - 9i{x{s), y{s))] ds 



\ 



[92{x{s), y{s)) - g2{x{s), y{s))] ds 



[gn{x{s),y{s)) - 5„(x(s), y(s))] ds j 
( T[gi{x{si),y{si)) - gi{x{si),y{si))] ^ 

T [52(2^(52), y(s2)) - 52(^(^2), y(s2))] 



\ T[5«(a;(s„),y(s„)) -g„(x(s„),2/(s„))] / 
Thus, using (|5.34p . one can obtain that 

[9{x(.s), y{s)) - g{x{s), y{s))] ds 

> r |gj„(x(sjj,2/(sjj) -5j„(x(sjo),2/(sjJ)| 
riieo 



> 



2n 



It is clear that, for t £ [9,9 + t], y{t) and y{t) satisfy the integral equations 



and 



y{t) = y{9)+ Ay{s)ds+ gix{s),yis))ds 



ft rt 
y{t)^y{9)+ I Ay{s)ds+ g{x{s),y{s))ds. 
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respectively, and herewith the equation 

y{t) - y{t) = {y{e) - y{e)) + f Mvi^) ~ y{s))ds 

Je 

+ [g{xis),y{s)) - g{x{s),y{s))]ds 
Je 

Je 

is achieved. Hence, we have the inequahty 



e+T 



\\y{e + T)-y{0 + T)\\ > 

-\\y{0)-ym 

^ \\A\\\\yis)-y{s)\\ds 

-T 

L3 \\y{s)~y{s)\\ ds. 



[gixis), y{s)) - g{x{s),y{s))]ds 



(5.35) 



Now, assume that max \\v(t) — v(t)\\ < ——r- tt— — n . mm ■ 

te[e,e+r]"^^ ^ - 2n[2 + t{L3 + \\A\\)] 



In the present case, one encounters with 



a contradiction since, by means of the inequahtics (|5.34[) and (j5.35p . we have 

max \\y{t)-y{t)\\ > hie + T)-y{e + T)\\ 

>^-[l + riL, + \\A\\)] max \\y{t) -ym 



tLiEo 1 + t{L3 + \\A 



2n[2 + r(L3 + ||A||)] 



2n \' 2 + t{L3 + \\A\\) 
rL leo 



2n[2 + r(L3 + P||)]- 
Therefore, one can see that the inequahty max \\y{t) — y{t)\\ > 



rLieo 



is vahd. 



te[e,8+r]"" ■ - "' 2n[2 + r(i3 + |jA||)] 
Suppose that at the point 77, the real valued function \\y{t) —y(t)\\ takes its maximum on the interval 
[9,9 + t], that is, 

max \\y{t) - y{t)\\ = \\y{r]) - y{7])\\ 
te[eM+T] 

for some 9 < 1] < 9 + t. 

For t E [9,9 + t], by favour of the integral equations 

y{t)^y{il)+j Ay{s)ds+ g{x{s),y{s))ds, 

J 7} J 7} 



and 



y{t) = yiv)+ Ay{s)ds+ g{x{s),y{s))ds, 

J 7) J T) 
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we obtain 



y{t) - vit) - iy{7]) - yir])) + / - yis))ds 



Define = 



+ / [9ixis),yis)) - gixis),yis))]ds. 



tL leo 



and lot fll = i - 
2'8n(MP||+i\/o)[2 + r(L3 + P||)]i ^ \ V - r\ if > 



that the interval = [6'^, 0^ + r^] is a subset of [0, 9 + t] and hence of J. 
For t e J^, we have 



? . We note 

2 



\yit)-yit)\\ > \\yiv)-yiv)\\ 



\A\\\\yis)-yis)\\ds 



\gixis),yis)) - gixis),yis))\\ ds 



V 



> 



> 



- 2M 



2n[2+T(L3+||A||)] 

2n^2+;k;||A||)] -2^^(^PII+Mo) 

rLicg 

4ri[2+T(L3+||A||)l ■ 



and the length of the interval 



Consequently, we get \\yit) — yit)\\ > ei, t £ J^, where ei 
does not depend on the functions a;(i),x(t) g ^4- 
The proof of the lemma is finaUzed. □ 

We say that the set is sensitive if both and s^y are sensitive. According to this definition, if si^x is 
sensitive, then Lemma |5 . 1 1 implies that the set ^ is also sensitive. 



4n[2+r(L3 + ||.4||)] 



Now, let us illustrate the replication of sensitivity through an example. It is known that the Lorenz 
system 



(T i — xi + X2) 
~X2 + rxi — Xix^ 
-~bx3 + X1X2, 



(5.36) 



with the coefficients a = 10, 5 = 8/3, r = 28 admits sensitivity [35]. We use system (|5.36p with the specified 
coefficients as the master system and constitute the 6— dimensional result-system 



10(-xi +a;2) 
—X2 + 28x1 — X1X3 

-|X3 + X1X2 
-5X4 + 2^3 

-2x5 + 0.0002(x2 - 
—3x6 — 3xi. 



(5.37) 



X5)3 +4X2 



When system (|5.37p is considered in the form of system (j2.10p + (|2.12l) . one can see that the diagonal 

matrix A whose entries on the diagonal are —5, —2, —3 satisfies that ||e"^* || < Ne~^*', where = 1 and lo = 2. 

p3 ^ 



We note that the function g : R'^ x R"* — > K'^ defined as 

5(xi, X2, X3, X4, X5, xe) = (x3, 0.0002(x2 - xs)^ + 4x2, 



-3xi 



provides the conditions (A4) and (A5) with constants Li = -^,-^2 = and L3 = | since the chaotic 

attractor of system (|5.37p is inside a compact region such that |x2| < 30 and |x5| < 50. Consequently, system 
(|5.37[) satisfies the condition (^7). 

In Figure El one can sec the 3— dimensional projections in the xi — X2 — X3 and X4 — X5 — xg spaces of 
two different trajectories of the result-system (j5.37p with adjacent initial conditions, such that one of them 
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is in blue color and the other in red color. For the trajectory with blue color, we make use of the initial 
dataa;i(0) = -8.57,X2(0) = -2.39, 2:3(0) = 33.08,a;4(0) = 5.32,X5(0) = 10.87,a;6(0) = -6.37 and for the one 
with red color, we use the initial data a:;i(0) = — 8.53,X2(0) = — 2.47,a;3(0) = 33.05, X4(0) = 5.33,a:;5(0) = 
10.86, a;6(0) = —6.36. In the simulation, the trajectories move on the time interval [0,3]. The results seen in 
Figure [9] supports our theoretical results indicated in Lemma 15.11 such that the replicator system, likewise 
the generator counterpart, admits the sensitivity feature. That is, the solutions of both the generator and 
the replicator given by blue and red colors diverge, even though they start and move close to each other in 
the first stage. 




Figure 9: Replication of sensitivity in the result-system 115. 371 1. (a) 3— dimensional projection in the xi — X2 — X3 space, (b) 
3— dimensional projection in the X4 — X5 — xg space. The sensitivity property is observable both in (a) and (6) such that the 
trajectories presented by blue and red colors move together in the first stage and then diverge. In other words, the sensitivity 
property of the generator system is mimicked by the replicator counterpart. 

Lemma 5.2 Transitivity oj s^x implies the same feature for s^y. 

Proof. Fix arbitrary small e > 0, arbitrary large E' > and let yify e s^y be a given function. Arising from 
the definition, introduced by ()2.15p . of the set s^y^ there exists a function x(t) e s^^ such that y(t) ~ 4'x{t){t)- 
Let 7 = 2MaN{L^-^NL^)+N L^Lu ' Condition {A7) guarantees that 7 is positive. Since there exists a solution 
x*{t) £ £/x, which is dense in one can find ^ > and an interval J C M with length E such that 
||a::(t) — x*{t + < 76, for all t E J. Without loss of generality, assume that J is a closed interval, that is, 
J ~ [a,a + E] for some real number a. 

Let y*{t) = <l>x*{t){t)- For t £ J, the bounded on M functions y{t) and y*{t) must satisfy the relations 

J —00 

and 

y*{t)^ f e^('-^^g{x*{s),y*{s))ds, 
J —00 

respectively. The second equation above implies that 

y*{t + 0= / e^^'+^-^^g{x*{s),y*{s))ds. 

J —00 
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t 

e" 



Using the transformation s = s — ^, and replacing s by s again, we get 

J —OO 

Therefore, for t G J, we have 

yit)-y*it + 0^ f e^(*-^)[5(x(s),2/(s))-5(x*(s + 0,2/*(s + 

J — OO 

s^(*-«) \g{x{s), y{s)) - g{x{s), y* {s + 0)]ds 
+ / e^(*-^)[5(x(s),y*(s))-5(x*(5 + e),y*(s + C))]ds, 

J a 

which imphes the inequahty 

\\y{t)-y*{t + 0\\< r 2MoNe-^^'-'Us 

J — OO 

+ /* iVLse-"^*-^) Ms) - y* {s + Oil ds 

J a 

+ f NL2e"^^'-'^\\x{s)-x*{s + ^)\\ds 

J a 

+ f NL^e"^^'-''^ \\y{s) ~y*{s + 0\\ ds. 

J a 

Hence, we have 

e"* hit) -y*it + Oil < ^-M^e- + (e^ - e^) 

UJ UJ 

+ / NL:,e^'\\y{s)-y*{s + 0\\ds. 

J a 

Through the implementation of Lemma 2.2 |9Q| to the last inequality, we attain 

e"* Mt) -y*it + Oil < ^M^I^e-^ + (g-* _ e^) 

'2MoN , NL2je 



UJ V ^ ^ J 



uj{uj - NL3) 
Multiplying both sides by e"'^*, we achieve that 

\\y(t)-y*it)\\ < ^M^eSNL3-^)(t-a) 



U! 

\ _ g(JVL3-")(*-a) 
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Now, suppose that E > In (^^"j . li t e [a+§,a + E], then it is true that e^^-^^-")^*-'^) < 

^{NL3-u:)§ ^ result, we have 



\\y{t)-y*it + o\\ < 



2MnN NLo 



UJ - NL3 
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for t <E Ji = [ai, ai + Ei] , where ai = a+-j and Ei = ^. Consequently, the set s^y is transitive in compliance 
with Definition 13.31 

The lemma is proved. □ 

Lemma 5.3 If s^x admits a dense countable collection of periodic functions, then the same is true for siy. 

Proof. Fix — (j>x(t){t) € s^y, arbitrary small e > and arbitrary large E > 0. Assume that ^4 
admits a dense countable collection of periodic functions. Let 7 = „» r nr^^^T/T'^N^i ^rr — • Note that 7 is 
a positive number by condition {AI). In this case, there exist x{t) £ and an interval J with length E 
such that \\x(t) — < 7e for t £ J. Without loss of generality, assume that J is a closed interval, that is, 
J = [a, a + i?] for some a G M. 

We note that by condition (^4) there is a one-to-one correspondence between the sets and '^y — 
{4'x(t){t) I x{t) £ 'Sx] , and xit) £ "^x is pT— periodic if and only if (f>x{t){t) € "^y is pT— periodic. Therefore 
'i^y C .s/y is a countable collection of periodic functions and our aim is to show that the set J^j, is dense in £/y 
in accordance with Definition 13.41 

Let y{t) ~ 0^(4) (t), which clearly belongs to the set '^^y. Making use of the relations 

t 



y{t)= e^^'-^^g{x{s),y{s))ds, 



and 



for i G J, we attain 



—00 



y{t)= e^^'-^'>g{x{s),y{s))ds, 



00 



y{t)~y{t)= / e^^'-^^[g{x{s),y{s))-g{x{s),y{s))]ds 

J —00 

e^^'-^^[g{x{s),y{s))~g{x{s),y{s))]ds 
' e^^'~^^[g{x{s),y{s))-g{x{s),y{s))]ds. 

'a 

The last equation implies that 

\\y{t)-ym < f 2Mo7Ve-"(*-^)ds 

J —00 

+ f NL^e-^('-''^\\y{s)-y{s)\\ds 



+ / iVLze""^*^'^ \\x{s) - x{s)\\ ds 



a; UJ ^ ' 

+ f NLse-'^^'-'^ \\y{s) - y{s)\\ ds. 
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Hence, we achive that 



\y{t)-ym < 



+ / NL^e'-'\\y{s)-y(s)\\ds. 



Application of Lemma 2.2 [90] to the last inequality provides us 



UJ \ 



+■ 



u}(lu — NL3) 

Multiplying both sides by e~'^*, we obtain 



NLst {u-NL3)t _ {0J-NL3) 



1 



ll!/(()-«()ll<— e'"'--"><'-"> 

UJ 



f NL2^e N'^L2L3^€ 



1 - e 



(JVL3-")(t-a) 



\ LU Uj{uj — NL3) ^ 



Suppose that E > In (^j^'j . U a + § < t < a + E, then one has e(^^3-'^)(t-a) < g(ArL3-^)f < ^g. 

'2MoN NL2 



U-NL3 

Consequently, the inequality 



\m-ym < 



OJ UJ — NL3 



76 = e, 



holds for ^ e Ji = [oi, oi + Ei] , where oi = a + -j and Ei = ^■ 
The proof of the Lemma is completed. □ 

Theorem 5.2 // the set ^4 is Devaney's chaotic, then the same is true for the set £/y. 



5.3 Machinery of Li-Yorke chaos extension 

In this part of our paper, the extension of the chaos in the Li— Yorke sense is taken into account. Our aim is 
to prove that if s/^ is chaotic in the sense of Li-Yorke, then the same is true for the set £/y and hence for the 
set £/. 

We start by presenting the following lemma which states the replication of proximality property. 
Lemma 5.4 If a couple of functions {x{t),x{t)) € s^x x is proximal, then the same is true for the couple 

{(t>x[t){t),4>x{t){t)) e X Sl/y. 

Proof. Fix arbitrary positive real numbers e and E. Let us define 7 = 2MoN^l^nl3)+n L2U) ^^ich is a 
positive number by condition {A7). Suppose that a given couple of functions {x{t),x{t)) g x is 
proximal. Therefore, one can find a sequence of real numbers {Ei\ satisfying Ei > E for each i S N, and 
a sequence {at} , — >■ 00 as i — >■ 00, such that we have \\x{t) —x{t)\\ < je, for each t from the intervals 

Ji = [ai,ai + Ei], i e N, and Ji C\ Jj = % whenever i ^ j. 
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Now, we shall continue our proof by fixing an arbitrary natural number i. Since the functions yit) 
4>x{t)it) € -^y a-nd y{t) = (j}x(t){t) G satisfy the relations 



y{t)^ I e^^'-^'>g{x{s),y{s))ds, 

and 

respectively, for t € Ji, we have 

yit) - y{t) - / e^(*-^) [g{x{s), y{s)) - g{x{s),y{s))]ds 



— OO 



+ / e^^'-^\g{x{s),y{s))~g{x{s),y{s))]ds 

J a-i 

+ f e^(*-^) [g{x{s)Ms)) - g{x{s),y{s))]ds. 

J a -i 

This implies that the inequahty 

\\y{t)-y{t)\\ < r 2AfoA^e-(*-^)ds 

^ — OC 

' NLse-'^^'~^'^\\y{s)-y{s)\\ds 



+ / NL2e-^<^'-'^\\x{s)-x{s)\\ds 



+ I NL:,e-"^'-'^\\y{s)-y{s)\\ds 

is valid. Hence, we achieve that 

e-*|ly(0-27Wll < 
+ f NL3e'^'\\y{s)^y{s)\\ds. 



OJ (jj 

t 



Implementing Lemma 2.2 I90| to the last inequality, we obtain 

e-* \\y{t) y{t)\\ < ^^e"- + (e^ - e"-) 



UJ OJ 



UJ V ^ ^ / 
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By multiplying both sides by e one can easily verify that 

2AIoN 



\y{t)-ym < 



^(NL3-u){t-a,) 



UJ uj{lj — NL:i) 



NL2je 



U! 



If E> 



u - NL-i 

In and t belongs to the interval [ai + ai 



(NL3-u)(t- 



Ei\ , then one has 



< 76. 

Since the natural number i was arbitrarily chosen, for each i £ N, we have 

'2MoN NL2 



for each t e J, 



ai,ai + 



\\y{t)~ym < 



where a. 



w - NL3 
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^ and Ei = Note that C K is an interval with length 



not less than i? = -j. As a consequence, the couple of functions {4>x{t)it)T <Px{t)it)) G s/y x si/y is proximal 



according to the Definition 13.61 

The proof of the Lemma is accomplished. □ 

Lemma 5.5 If a couple of functions {x{t),x{t)) e x •s/x is frequently (eg, A) —separated for some positive 
real numbers eo and A, then the couple of functions {(j)x{t) (^)j 'fe(t)(^)) ^ £^yX£/y is frequently (ei, A) — separated 
for some positive real numbers ei and A. 

Proof. Suppose that a given couple of functions {x{t),x{t)) e s^x x ^4 is frequently (eo, A) separated, for 
some eo > and A > 0. In this case, there exist infinitely many disjoint intervals, each with length not less 
than A, such that \\x{t) — x{t)\\ > eq, for each t from these intervals. Without loss of generality, assume that 
these intervals are all closed subsets of R. In that case, one can find a sequence {A^} satisfying A^ > A, 
i e N, and a sequence {di} , — > 00 as i — > 00, such that for each i £ N the inequality \\x{t) — x{t)\\ > eo 
holds for t £ Ji ^ [di,di + Ai], and Ji n Jj = whenever i ^ j. 

Throughout the proof, let us denote y{t) = (j)x{t){t) £ s^y and y{t) = (/%(t)(i) £ -s^y 

Our aim is to show the existence of positive real numbers ei,A and infinitely many disjoint intervals 
Ji C Ji,i £ N, each with length A, such that the inequality \\y{t) —y{t)\\ > ei holds for each t from the 
intervals Ji,i £ N. 

( 9i{x,y) \ 



As in the proof of Lemma l5. 11 we again suppose that g{x,y) 



92ix,y) 



V gnix,y) J 



where each gj, 1 < j < n, 



is a real valued function. 

Using the equicontinuity on K of the family defined by equation (|5.3ip . one can find a positive real 
number r < A, independent of the functions x{t),x{t) £ £/x, y{t),y{t) £ sz/y, such that for any ti,i2 £ K 
with |ti — < T the inequality 



1(5, {x{t,),y{t,))-gj (x(ii),2/(ii))) - {g, (.T(t2), 2/(^2)) - 5, (x(t2), 2/(^2)))! < 



2n 



(5.38) 



holds, for all 1 < j < n. 
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Suppose that the sequence {s;} denote the midpoints of the intervals J^, that is. Si = di + ^ , for each 
i G N. Let us define a sequence {6i} through the equation Oi ^ Si — ^. 
Now, let us fix an arbitrary natural number i. 

In a similar way to the method specified in the proof of Lemma 15.11 one can show the existence of an 
integer ji — ji{si), 1 < ji < n, such that 



On the other hand, making use of the inequality ()5.38|) . it is easy to verify that 

\9]^ {x{si),y{si)) ~ gj, (x(s,),y(s,))| - \gj^ {x{t),y{t)) - g^^ {x{t),y{t))\ 

< \i9j^ - 9j. {xit),y{t))) - {gj, (x(s,),y(s,)) - gj, (x(s,), y(s,)))l 



(5.39) 



< 



2n ' 



for all t £ [Oi,di + t] and therefore by favour of (|5.39p . we obtain that the inequality 

Ife ixit),yit)) - gj^ (x{t),y{t))\ > \gj^ {x{s,),y{si)) - g,^ (x(si), y(s,))| - > 



2n 



2n 



(5.40) 



is valid on the same interval. 

Using the mean value theorem for integrals, it is possible to find real numbers s\, s^, . . . , sl^ & [9i, 9i + t] 
such that 



h+r 



[g{x{s), y{s)) - g{x{s), y{s))] ds 



[gi{x{s), y{s)) - gi{x{s), y{s))] ds 



[92{x{s), y{s)) - g2{x{s), y{s))] ds 



[g„(x(s), y(s)) - gn{x{s), y{s))] ds 



'gi{x{s\),y{s\)) - gi{x{s\),y{s\)) 
92{x{s\),y{s\)) - g2{x{s\), y{s\)) 



\ 



\ [5"(a;(4),y(4)) -5„(x«),y«))] / 
Hence, the inequality ()5.40p gives us that 

[9{x(.s), y{s)) - g{x{s),y{s))] ds 



> T\gj^{x{sjJ,y{s]J) - gj^{x{s)J,yis].))\ 
tLiCo 



> 



2n 



For t € [6i, 9i + t], the functions y{t) G s^y and y{t) G s^y satisfy the relations 



y{t)=y{e,)+ Ay{s)ds+ g{xis),yis))ds 
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and 



f* ft 
y{t)^y{e,)+ I Ay{s)ds+ g{x{s),y{s))ds, 



respectively, and herewith the equation 



y{t) - y{t) = {y{e.,) - y{e.,)) + / A{y{s) ~ y{s))d. 



+ [9ix{s),y{s)) - g{x{s),y{s))]ds 



[gi^is), y{s)) ~ g{x{s), y{s))]ds 



is achieved. Taking i = + r in the last equation, we attain the inequality 



\M+r)-y{9, + T)\\ > 



[g{x{s), y{s)) - g{x{s), y{s))]ds 



\m)-ym\- {\\A\\+L,)\\y{s)-y{s)\\ds 



(5.41) 



Now, assume that max \\v(t) — y(t)\\ < ——r- tt— — ttttt- 

tion since, by means of the inequalities (j5.40p and (j5.4ip . we have 



In this case, one arrives at a contradic- 



max \\yit)-yit)\\ > l|y(^. + ^) - ^ + ^)ll 
te[9i,ei+T] 

>^-[l + riL, + \\A\\)] max Mt) - yit)\\ 



>^-[l + r{Ls + \\A\\)] 



tLiCo 



2n[2 + t{L3 + \\A\\)] 

= ^-^1^0 (l _ 1 + t(£3 + P||) \ 
2n \ 2 + TiLs + \\A\\) J 

rLiep 

" 2n[2 + T(L3 + PI|)]' 

Therefore, it is true that max lluft) — > — ; 7—— — ,, , i,,, . 

tG[e.,s.+T]"^^ ^ '^^ ^" 2n[2 + r(L3 + ||A||)] 

Suppose that the real valued function ||y(t) — assumes its maximum on the interval [0^, Qi + t], at a 
point r\i. In other words, for some r\i £ [6'i, 6*^ + t] , we have 



max \\y{t) - y[t)\\ = \\y{m) - y{m)\\ 



Making use of the integral equations 



yit)=yiVi)+ / Ay{s)ds+ / g{x{s),y{s))ds, 



and 



y{t)=y{Vi)+ Ay{s)ds+ g{x{s),y{s))ds, 

J j-ji J rii 
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on the time interval [6i , 6i + t] , one can obtain that 



Define A = min 



y{t) - y{t) = {y{n,) - y{i^,)) + / A(y(s) - y{s))ds 

Jvi 

t 

\g[x{s),y{s)) - g{x{s),y{s))]ds. 

tLiCq 



2 ' 8n(M \\A\\ + Mo)[2 + t(L3 + 1|A||)] 



and let 9} = 



if 77 < 9, 
if Vt > 



For each t € [9], 9} + A], we have 

\\y{t)-ym>Mv^)-y{v^)\\- 



\A\\\\y{s)-y{s)\\ds 



> 



> 



g{x{s), y(s)) - g{x{s), y{s))\\ ds 
2M \\A\\ ri - 2A/ot1 



2n[2+r{L3 + lt^l|)] 

2n[2H-;j;W||)] -2r^(^^PII+M^ 
4nl2+r(L3 + l|A||)l ' 



The information mentioned above is true for arbitrarily chosen the natural number i. Therefore, for each 
i E N, the interval J; = [91,9} + A] is a subset of [9i,9i + r], and hence of Ji. Moreover, for any i S N, we 
have \\y{t)-y{t)\\ >ei,te where ei ^ 4„[2+r(23+||A||)] • 

Consequently, according to Definition 13. 7[ the couple of functions {4>x{t){i)i4>x{t){'t')) G 'S^y x ■^y is fre- 
quently (ei, A)— separated. 

The proof of the lemma is finalized. □ 

We end up the present subsection by stating the following theorem and its corollary. 
Theorem 5.3 // the set s^x is Li-Yorke chaotic, then the same is true for the set s^y. 



Proof. Assume that the set s^x is Li-Yorkc chaotic. As we indicated in Subsection l5.1[ it is easy to show that 
for any natural number k, x(t) G is a fcT— periodic function if and only if 4>x{t) (t) G is also fcT— periodic, 
where "i^x and denote the sets of all periodic functions inside jz/x and £/y, respectively. Therefore, the set 
£/y admits a /cT— periodic function for any fc € N. 

Now, suppose that the set '^x is a scrambled set inside £/x and define the set 



% ^ {(l>xit){t) \x{t)e'^x}- 



(5.42) 



Condition (^4) implies that there is one-to-one correspondence between the sets 't^x and 'T^y. Since the 
scrambled set "t^x is uncountable, it is clear that the set "^j, is also uncountable. Moreover, using the same 
condition one can show that no periodic functions exist inside since no such functions exist inside "^x- 



That is, % n % 



Since each couple of functions inside "^ipx x '^x is proximal, Lemma 
couple of functions inside '^y x . 

Similarly, Lemma [^3] implies that if each couple of functions {x{t),x{t)) g x 



implies the same feature for each 



fx) is frequently 



(eo. A)— separated for some positive real numbers eo and A, then each couple of functions {y{t),y{t)) G'^yX 
C^y X ^y) is frequently (ei. A)— separated for some positive real numbers ei and A. Consequently, the set 
is a scrambled set inside £/x, and according to Definition 13.101 s/y is Li-Yorke chaotic. 

The proof of the theorem is accomplished. □ 

An immediate corollary of Theorem 15.31 is the following. 



Corollary 5.2 If the set six is Li-Yorke chaotic, then the set si is chaotic in the same way. 
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5.4 Consecutive replications of chaos 

In this part of the paper, in addition to the system (|2.10p + (|2.12p . we take into account the system 

z' Bz + h{y,z), (5.43) 

where h : M" x K' K' is a continuous function in aU of its arguments, and the constant I x I real valued 
matrix B has real parts of eigenvalues all negative. 

It is easy to verify that there exist positive real numbers TV and uj such that ||e^*|| < iVe~"*, for all t>0. 

In our next theoretical discussions, the system (|5.43p will serve as the second replicator system in the 
mechanism presented by Figure [31 and we need the following assumptions which are counterparts of the 
conditions (A4) — {A7) presented in Section [2j 

(A8) There exists a positive real number Li such that \\h{yi, z) — h{y2, z)\\ > Li \\yi — y2\\ , for all yi, y2 G M", 
z e R'; 

(A9) There exist positive real numbers L2 and L3 such that \\h(yi,z) — h{y2,z)\\ < L2 \\yi — y2\\ , for all 
yi,y2 e M", 2 e M', and \\h{y,zi)-h{y,Z2)\\ < ZgHzi -Z2II , for aU y e R", zi,Z2 G R'; 

(AlO) There exist a positive number Kq < 00 such that sup \\h{y, z)\\ = Kq] 
(All) NL3-U < 0. 

Likewise the definition for the set of functions £/y, given by (|2.15p . let us denote by jz/z the set of all 
bounded on K solutions of system z' = Bz + h{y{t),z), for any y{t) £ s^y. 
In a similar way to the Lemma p.lj) . one can show that the set 

= {z(t) I z(t) is a solution of the system z' =- Az + g{y{t), z) for some y{t) G '^y } . (5.44) 

is a basin of £/z ■ 

Similar results of Theorem 14.11 introduced in Section |4l and of the Theorems 15. li 15.21 and 15.31 presented 
in Subsections 15. li 15.21 and 15. 3i respectively, hold also for the set ia4 , and we state this result in the next 
theorem. 

We note that, in the presence of arbitrary finite number of consecutive replicator systems, which obey 
conditions that are counterparts of (^4) — {A7), one can prove that a similar result of the next theorem also 
holds. 

Theorem 5.4 // the set si^ is chaotic through period- doubling cascade or Devaney chaotic or Li-Yorke 
chaotic, then the set ^/^ is chaotic in the same way as both ^ and £^y. 

Proof. In the proof, we will show that for each z[£) G s^z and arbitrary (5 > 0, there exist z(t) G s^z and 
to G M such that ||z(to) — ^(io)|| < (5, which is needed to show sensitivity of The remaining parts of the 
proof can be performed in a similar way to the proofs presented in Subsections 15 . li W72\ and and therefore 
are omitted. 

Suppose that the set si^ is sensitive. Fix an arbitrary (5 > and let z(t) £ £/z he a. given solution of 
system (|5.43p . In this case, there exists y{t) = 4'x{t}{t) G £/y, where x{t) G £/xi such that z{t) is the unique 
bounded on M solution of the system 2' = Bz + h{y{t), z). 

Let us choose a positive number e — e{6) small enough such that the inequality 

/ NL2 \ f NL2 \ 
1 + ^ 1 + 771- e< (5 

\ 2} - A^ia ) \ ^- NL3 J 

holds, and denote ei = (^1 + lj^nLs ) ^' ^^''^i take R = R{e) < sufficiently large in absolute value 
such that both the inequalities MkKe-iNL,-u;)§ < ^ .^^^ M4iiVg-(A'L3-S)f < ^re valid, and let 61 = 
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Si{e,R) — ee^"^ . Since the set ^ is sensitive, one can find x{t) G jz^ and to G M such that the inequality 
\\x{to) -x{to)\\ < Si holds. 

As shown in the proof of Lemma [5Jl for t e [to + R, to], one has 

\\xit)-xm <e, 

and 

LU ~ I\ 1^3 UJ 

As a consequence of the last inequality, we have 

\\yit)-y{t)\\<€i, 

for te[to + §,to]. 

Suppose that z(t) is the unique bounded on M solution of the system z' = Bz + h{y{t),z). One can see 
that the relations 

z{t) = f e^('-'^h{y{s),z{s))ds 



and 



are valid. 
Therefore, 



z{t)= f e^^'-^h{y{s),z(s))ds, 



z{t) - z{t) = / e^(*-^) [h{y{s), z{s))^ h{y{s),z{s))]ds 

J — OO 

and hence we obtain 

II^W - z{t)\\ < f ^e-"(*-^) \\h{y{s), z{s)) - h{y{s) ,z{s))\\ ds 
+ f Ne-~-('-^^ \\h{y{s),z{s)) - h{y{s),^{s))\\ ds 
iVe--(*-^) \\h{y{s), z{s)) - z(s))|| ds. 

-OO 

Since \\y{t) ~y{t)\\ < t\ for t G [to + -^,^0], one has 



it) - z(t)\\ < NL3 / e-=(*-^) ||z(s) - z{s)\\ ds 

Jto + § 
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Now, introduce the funetions u{t) ~ e"* ||z(i) — z(t)|| , k{t) 

( 2MoN _ JVZaei A C(to + f ) 
I a; ^ J 

These definitions give us the inequahty 



NL2ei „ut 



J-e"*, and v{t) = c + k{t) where 



u{t) < v{t) + / NL3u{s)ds. 

Jto + § 

Applying Lemma 2.2 [90] to the last inequality, one can verify that 



u{t) < v{t) + NL3 / e^-^='(*~'*)/i(s)ds. 



to+-. 



Therefore, for t £ [to + ^,ta] we have 



Hence, 



< c + kit) + c (e^£3(t-*o-f ) _ 1) + N^Lg^^^z.t 



kit) + ce^23(t-*o-f ) + jV^^2£3^i 

— NL^) 



1 „ g(^Z3-2)(t-to--|) 



1 - e 



(AfL3-S)(t-to-#) 



,(ArL3-2)(t-to-f ) 



Consequently, 



+ 



N^L2L3ei 
ojiuj — NL3) 



w - NL3 

NL2ei 
uj - NL^ 



I _ g(A'L3-C)(t-to-f ) 



1 _ g(wZ3-2)(t-to-f ) 



i^l^g(ArL3-2)(t-to-f )^ 



\\zito)~zito)\\ < %^ + ^^^(^-NL,)- 



UJ - NL3 



< 1 



<S. 



ei 



The theorem is proved. □ 



6 Controlling morphogenesis of chaos 

In previous section we have theoretically proved self-replication of chaos for specific types. Controlling the 
self-replicated chaos is another interesting problem and in this section we give an answer. The next theorem 
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and its corollary indicates a method to control the chaos of the replicator system (|2.12|) and the result-system 
(|2.10p + (|2.12p . respectively, and reveals that controlling the chaos of system (|2.10p is sufficient for this. 

Theorem 6.1 Assume that for arbitrary e > 0, a periodic solution Xp(t) S £/x "is stabilized such that for any 
solution x{t) of system (|2.10p there exist real numbers a and E > such that the inequality \\x{t) — Xp{t)\\ < e 
holds for i e [a,a + E]. 

Then, the periodic solution 4>xp{t) (t) € -s/y is stabilized such that for any solution y{t) of system (j2.12p there 

exists a number b > a such that the inequality \\y{t) — '/'a:p(t)(i)|| (-'^ ^ uj'^nLs ) ^ holds for i G [6,a + E], 
provided that E is sufficiently large. 



Proof. Fix an arbitrary solution y{t) of system y' = Ay + g{x{t),y) for some solution x{t) of system (|2.10p . 
According to our assumption, there exist numbers a and E > such that the inequality ||a;(i) — a:;p(t)|| < e 
holds for t £ [a, a + E]. Let us denote yp{t) = (j)^ (t)(0 € •^y It is clear that the function yp{t) is periodic 
with the same period as Xp(t). Since y{t) and yp{t) satisfy the integral equations 



and 



respectively, one has 



yit) = e^(*-'^)y(a) 
yp{t) = e^^'-^^yp{a) 



s^(*-^),9(.T(s),y(s))ds, 



9{xpis),ypis))ds, 



y{t)^yp{t)^e^^'-^\y{a)-yp{a)) 



[9{x{s),y{s)) - 9{xp{s),yp{s))] ds 



\y{a) - yp{a)) 
,^^'-'^[g{x{s)M)-9{<s),yp{s))]ds 

[g{x{s),yp{s)) - g{xp{s),yp{s))] ds. 
By the help of the last equation, we have 



,A(t-a) 



\\y{a)-yp{a) 











g{x{s),yp{s))\\ ds 


[ 


e^(*-^) 






- ff(a;p(s),yp(s))|l ds 



\\y{t)-yvm < 
+ 

+ 

<^re-(*-)||y(a)-y,(a)|| 
+ f Ne--^'-^^L,\\yis)~yp{s)\\ds 

Ne-'^^'-''>L2\\x{s)- Xp{s)\\ds 
< Ne-^'-^U\y{a) - yp{a)\\ 
) 

+NL3 I e-^^'-''>\\y{s)-yp{s)\\ds. 
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Let u : [a,a + E] ^ R+ be a function defined as u{t) ~ e"* \\y{t) — yp{t)\\ . Then we reach the inequality 



) +NL3 / u{s)ds. 



Now, define ip : [a,a + E] ^ R+ as tp{t) = EL^i^'^t and similarly cj) : [a,a + E] ^ R+ as (j){t) = + c 
where c = Ne^'^ ^ J/pl*^)!! ^ EL2l^<^°-_ gy means of these definitions we get 



u{t) < <j){t) + 7VL3 / u(s)ds. 

a 

Implementation of Lemma 2.2 [90| to the last inequality, where < G [a, a + provides us the inequality 
u{t) <c + iit) + NL3 f e^^'^'-'^cds + NLs [ e^-^^(*-"V(s)ds 



and therefore, 



U{t) <C + 1p{t) + C (^e'^Ls{t-a) - 1^ + 

NL2e 



lj{lu — NL^) 
7V||y(a)-2;p(a)||e""e^^^(*-'^) 



^1 _ ^{NL3-u:)(t~-a)\^ 



Consequently, 



< 



NL2e 



+ iV||2/(a)-2;p(a)|le(^^-^-)(*-'') 



NL2e 



— NL3) 



iV||2/(a)-2;p(a)||e(^^-^-)(*-'^) 



NL2e 



1 - e 



{NL3-ui){t-a) 



u) — N L3 \ / 

<iV||2/(a)-2;p(a)||e(^^^-")(*-'^) + 



a; - NL3 
NL2e 

LO - NL3 ' 



If y{a) = yp{a), then clearly \\yp{t) - y{t)\\ < (l + ^^nl^ ) ^ ^ [a,a + £']. Suppose that y(a) 7^ 
yp{a). For t > a + ^^1^^ In ( ^n,(.)l,^(.)i| ) , the inequality e(^^3-^)(*-a) < 
sume that the number E is sufficiently large so that E > 



N\\y{a)-yp{a) 



is satisfied. As- 



{ 



max < a.a + 



■In 



{a)-Vp{a)\\ ) } 



Thus, taking 6 



and i? = min jfi", i? — ^^r"^ hi 



WLa-i.^ I JV||y(a)-ap( 



a)ll)} 



one attains that 



ui-NLr 



e, for t E [b,b + E]. Here the number E stands for the duration of control for 



system (|2l^ . We note that b > a, < E < E and b + E = a + E. 



NL2 



e, for t e[b,a + E]. 



Hence - yp(t)|| < (^1 , ^^^^^ 
The proof is finalized. □ 

An immediate corollary of Theorem 16.11 is the following. 
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Corollary 6.1 Assume that the conditions of Theorem K.ll hold. In this case, the periodic solution Zp{t) = 
(a;p(i), 03.^(4-) (t)) G £/ is stabilized such that for any solution z(t) of system (|2.10p + (j2.12p there exists a 

number b > a such that the inequality \\zp{t) — z{t)\\ < ^2 + ^ e holds for t € [b,a + E], provided that 

E is sufficiently large. 

Proof. On account of the inequality \\z{t) ~ Zp{t)\\ < \\x{t) — Xp{t)\\ + \\y{t) — (l>xp{t)it)\\ i and using the 
conclusion of Theorem 16.11 one can show that \\zp(t) - z{t)\\ < (2+ ^ e holds for t G [b,a + E] and for 

some b > a. 

Remark 6.1 As a conclusion of the Theorem ] 6. 11 the transient time for control to take effect may increase 
and the duration of control may decrease as the number of consecutive replicator systems increase. 

In the remaining part of this section, our aim is to present an illustration which confirms the results of 
Theorem l6.1l and for our purposes, we will make use of the Pyragas control method |91] . Therefore, primarily, 
we continue with a brief explanation of this method. 

A delayed feedback control method for the stabilization of unstable periodic orbits of a chaotic system 
was proposed by Pyragas [91j . In this method, one considers a system of the form 

x' = H{x,q), (6.45) 

where q = q{t) is an externally controllable parameter and for (7 = it is assumed that the system (|6.45p is in 
the chaotic state of interest, whose periodic orbits are to be stabihzed [46l [87l [911 192] . According to Pyragas 
method, an unstable periodic solution with period ^ of the system (|6.45p with q = 0, can be stabilized by 
the control law q{t) = C [s (i — ^) — s{t)] , where the parameter C represents the strength of the perturbation 
and s(t) ~ a [x{t)] is a scalar signal given by some function of the state of the system. 

It is indicated in that in order to apply the Pyragas control method to the chaotic Duffing oscillator 

a;" + 0.10a;' - {l - x"^) = 0.24sin(t), (6.46) 

one can make use of the new variables xi = x,X2 = x' ,X3 = t, to construct the system 

X[ = X2 

x'2 = -0.10x2 + ixi (1 - xj) + 0.24sin(a;3) + C [x2{t - ^) - X2{t)] (6.47) 
4 = 1, 

where q{t) = C [x2{t — £,) — X2{t)] is the control law and the less unstable 27r— periodic solution can be 
stabilized by choosing the parameter values C = 0.36 and ^ = 27r. 

Now, we shall indicate an implementation of Theorem 16.11 through simulations. Let us combine system 
()6.47|) with two consecutive replicator systems and set up the 9— dimensional result-system 

X'l = X2 

4 = -0.10x2 + (1 - xl) + 0.24sin(a;3) + C [x2{t - - X2it)] 
4 = 1 

X4 = — O.lxi 

= -3x4 - 2a;5 - 0.008x1 + 1.6x2 (6.48) 

x'e = -X6 + X3 

X7 = X8 + 0.6x4 

x's = -3.1x7 - 2.1x8 - O.OOTxf + 2.5x5 

Xg = — Xg + 0.2X6. 

Since our procedure of morphogenesis is valid for specific types of chaos such as in Devaney's and Li- 
Yorke sense and through period-doubling cascade, we expect that our procedure is also applicable to any 
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other chaotic system with an unspecified type of chaos. Accordingly, the system (|6.48p is chaotic for the 
parameter value C = 0, since the system (|6.47l) admits the chaos for the same value of the parameter. 

Let us consider a solution of system (|6.48p with the initial data a;i(0) = X5(0) = xs{0) = 0.1,a;2(0) = 
-0.8,2:3(0) 0,2:4(0) = -0.5,2-6(0) = -0.1,2:7(0) = -0.2,2:9(0) = -0.3. We let the system evolve freely 
taking C = until t ~ 70, and at that moment we switch on the control by taking C = 0.36. At time t = 200, 
we switch off the control and start to use the value of the parameter C = again. In Figure [10] one can 
see the graphs of the 2:2, 2:5, coordinates of the solution. After switching off the control, the 27r— periodic 
solution loses its stability and chaos emerges again. 
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Figure 10: Pyragas control method applied to the result-system lj6.48|l . (a) The graph of a;2— coordinate, (b) The graph of 
0:5— coordinate, (c) The graph of xg— coordinate. The result of Pyragas control method applied to the generator system H6.47|l 
is seen in (a). Through this method, the 27r— periodic solution of the generator and accordingly the 27r— periodic solutions of 
the first and the second replicator systems are stabilized. In other words, the morphogenesis of chaos is controlled. For the 
coordinates xi,Xi and we have similar results which are not just pictured here. The control starts at f = 70 and ends at 
t = 200, after which emergence of the chaos is observable again. 



7 Discussion 

In this part of our paper, we intend to consider not rigorously proved, but interesting phenomena which can 
be considered in the framework morphogenesis of chaos. So we shall discuss again obtained above results and 
say about the possibility of the morphogenesis for other type of chaos, for intcrmittcncy, for relay systems, 
and discuss the possibility of quasiperiodic motions as an infinite basis of chaos. 

The type of chaos for the double-scroll Chua circuit is proposed in paper [93j . It is an interesting problem 
to prove that this type of chaos can be replicated through the method discussed in our paper. Nevertheless, 
let us show by simulations that the regular behavior in replicator systems can be also seen. This means that 
next special investigation has to be done. 

Moreover, this will show how one can use morphogenesis not only for chaos, but also for Chua circuits by 
uniting them in complexes in electrical (physical) sense, and observing same properties as unique separated 
Chua circuit admits. This is an interesting problem which can give a light for the complex behavior of huge 
electrical circuits. 

We start our discussions with morphogenesis of intermittency. 
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7.1 Morphogenesis of intermittency 

In Section [5l we have rigorously proved self-replication of specific types of chaos such as period-doubling 
cascade, Devaney's and Li-Yorke chaos. Consequently, one can expect that the same procedure also works 
for the intermittency route. 

Pomeau and Manneville j62] observed chaos through intermittency in the Lorenz system (|5.36p , with the 
coefficients a = 10, b = 8/3 and values of r slightly larger than the critical value ~ 166.06. To observe 
intermittent behavior in the Lorenz system, let us consider a solution of system (j5.36p together with the 
coefficients cr = 10, 6 = 8/3, r = 166.25 using the initial data xi(0) = —23.3, 2:2(0) ~ 38.3 and 13(0) — 193.4. 
The trajectories of the xi, X2 and 2:3 coordinates of the solution arc indicated in Figure [TTl where one can 
see that regular oscillations are interrupted by irregular ones. 




Figure 11: Intermittency in tlie Lorenz system 115.3611 . wliere cr = 10, b = 8/3 and r = 166.25. (a) Tlie graph of tlic 
xi— coordinate, (b) The graph of the a;2— coordinate, (c) The graph of the X3— coordinate. 

To perform morphogenesis of intermittency, let us consider the Lorenz system (|5.36p as a generator and 
set up the 6— dimensional result-system 

x[ = a (-xi + X2) 

x'2 ~ -~X2 + rxi — a;ia;3 

4 = -6X3 + XlX2 ,^ 

x'/^ = —X4 + Axi 

x'^= X(, + 2X2 

x'q = -3x5 - 2xe - O.OOOOSxi + 0.5x4, 

again with the coefficients cr = 10, & = 8/3 and r = 166.25. It can be easily verified that condition {A7) is valid 
for system (|7.49|) . We consider the trajectory of system (|7.49p corresponding to the initial data xi (0) = —23.3, 
2:2(0) = 38.3, X3(0) = 193.4, 0:4(0) = -17.7, 2:5(0) = 11.4, and a:6(0) ^ 2.5, and represent the graphs for 
the 0:4, 2:5 and xq coordinates in Figure [T^ such that the intermittent behavior in the replicator system is 
observable. The similarity between the time-series of the generator and the replicator counterpart shows the 
morphogenesis of intermittency. 

In next our discussion, we will emphasize by means of simulations the morphogenesis of the double-scroll 
Chua's attractor in a unidirectionally coupled open chain of Chua circuits. Approaches for the generation of 
hyperchaotic systems have already been discussed making use of Chua circuits which are all chaotic [TU |M] . 
It deserves to remark that to create hyperchaotic attractors in previous papers, others consider both involved 
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Fi gur6 12: Intermittency in the replicator system, (a) The graph of the 0:4— coordinate, (b) The graph of the X5 — coordinate, 
(c) The graph of the ig— coordinate. The analogy between the time-series of the generator and the replicator systems indicates 
the morphogenesis of intermittency. 

interacting systems chaotic, but in our case only the generator is chaotic and aU other Chua systems considered 
as repUcators are non-chaotic. 



7.2 Morphogenesis of the double-scroll Chua's attractor 

There is well known result of the chaoticity based on the double-scroll Chua's attractor [53] . It was proven 
first in the paper [53] rigorously and the proof is based on the Shil'nikov's theorem Since the Chua 

circuit and its chaotic behavior is of extreme importance from the theoretical point of view and its usage 
area in electrical circuits by radio physicists and nonlinear scientists from other disciplines, one can suppose 
that morphogenesis of the chaos will also be of a practical and a theoretical interest. 

We just take into account a simulation result which supports that morphogenesis idea can be developed 
also from this point of view. 

Let us consider the dimensionlcss form of the Chua's oscillator given by the system 



x[ = ka[x2 - xi ~ f{xi)] 
x'2 = k{xi — X2+ xz) 
= k{~(ix2 -~ 72:3) 
fix) = hx + Q.b{a-h) {\x ■ 



(7.50) 



II 



1|) 



21.32 
5.75 



, /? = 7.8351, 7 



1.38166392 



where a, 13,^, a,b and k are constants. 

In paper [97] , it is indicated that system (|7.50|) with the coefhcicnts a 
a = —1.8459, b ~ —0.86604 and k = 1 admits a stable equilibrium. 

In what follows, as the generator, we make use of the system (|7.50p together with the coefficients a = 
15.6, (3 = 25.58, 7 = 0, a = ~f , ^ = and k — 1 such that a double-scroll Chua's attractor takes place [98] 
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and consider the 12— dimensional result-system 





15.6[a;2 - fxi + ^ (jxi + 1| + |xi - 1|)] 








A = 


Xi~ X2+ X^ 








4 = 


-25.58x2 










2ji-3g2 0.13396.T4 + 0.48993 (|x4 + 1| 


+ \Xi 




+ 2x1 


4 = 


Xi — X^ + X(, + 5X2 








4 = 


7.8351X5 1-38166392 _^ 2X3 








4 = 


2ji-3g2 0.13396x7 + 0.48993 (|x7 + 1| 


+ \X7 




+ 2X4 


Xg = 


Xy — Xg + Xg + 3X5 








Xg = 


7.8351X8 1-38166392 Q.OOlXg 








^10 ~ 


= 2^1-32 [xn 0.13396x10 + 0.48993 (|xio ^ 




|a;io - 


I|)]+4x7 


^11 = 


= Xio — Xu + Xi2 - 0.1X8 








^12 - 


-- 7.8351X11 1-38166392^^^ _^_22,g. 









(7.51) 



In Figure [131 we simulate the 3— dimensional projections on the Xi — X2 — X3 and X4 — X5 — Xg spaces 
of the trajectory of the result-system (|7.51[) with the initial data Xi(0) = 0.634, X2(0) = —0.093, X3(0) = 
-0.921, X4(0) ^ -8.013, X5(0) = 0.221, X6(0) = 6.239, X7(0) = -50.044, X8(0) = -0.984, X9(0) = 48.513, 
a;io(0) = —256.325, xii(O) = 7.837, xi2(0) = 264.331. The projection on the xi — X2 — X3 space sho-ws the 
double-scroll Chua's attractor produced by the generator system (|7.50p . and projection on the X4 — X5 — Xg 
space represents the chaotic attractor of the first replicator. 




Figure 13: 3— dimensional projections of the chaotic attractor of the result-system 117.511 1. (a) Projection on the xi — X2 — 
space, (b) Projection on the x^ — x^— xq space. The picture in (a) shows the attractor of the original prior chaos of the generator 
system 1 17. SDH and (6) represents the attractor of the first replicator. The resemblance between shapes of the attractors of the 
generator and the replicator systems makes the morphogenesis of chaos apparent. 

In a similar -way, -we display the projections of the same trajectory on the X7 — X8 — xg and xio — xn — X12 
spaces, -which correspond to the attractors of the second and the third replicator systems, in Figure (Mj and 
the projections on the X3 — xg — xg and xi — X5 — xio spaces in Figure [TSl The illustrations sho-wn in Figure 
[131 Figure [HI and Figure [TS] indicate the morphogenesis double-scroll Chua attractor. 

Possibly the result-system (|7.5ip produces a double-scroll Chua's attractor -with hyperchaos -where the 
number of positive Lyapunov exponents are more than one and even four. 
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Fi gUre 14: 3 — dimensional projections of the chaotic attractor of the result-system (17.51I I. (a) Projection on the — xs — xg 
space, (b) Projection on the xio — xn — X12 space, (a) and (6) demonstrates the attractors generated by the second and the 
third replicator systems, respectively. 



7.3 Morphogenesis and the logistic map 

Since the logistic map is very convenient to verify Devaney's, Li-Yorke chaos and also chaos through period- 
doubling cascade, one can expect that it can be used as a generator. This is true in the theoretical sense, but 
not for simulations since only chaos through period-doubling route can be presented through simulations. 
In paper |44] it is rigorously proved that the Duffing equation perturbed with a pulse function 

x" + O.lSx' + 2x + 0.00004x3 = 0.02 cos(27rt) + to, /ioo), (7.52) 

with the coefficients mo = 2, mi = 1 and /ioo = 3.8, where the relay function iy{t,tQ, fioo) is described by 
equation (|5.22p . admits the chaos through period-doubling cascade on the time interval [0, 00) and obeys the 
Feigenbaum universal behavior [5S] . 

By favour of the new variables xi = x and X2 = x' , equation (|7.52p can be reduced to the system 

x'^ = -0.18x2 - 2x1 - 0.00004x? + 0.02cos(27ri) + i/(i,to,Aioo). ^ ' 

One can see that system (|7.53p is in the form of system (|5.21|) . For the illustration of morphogenesis of chaos, 
we will make use of the relay-system (|7.53p as the generator, in the role of the core as displayed in Figure 21 
and attach three replicator systems to obtain the 8— dimensional result-relay-system 

x'l = X2 

x'2 = -0.18x2 - 2x1 - 0.00004xf + 0.02 cos(27rt) -f v{t, to, Moo) 
X3 = X4 — O.lxi 

X4 = -10x3 - 6x4 - O.OSxj^ -I- 4x2 



(7.54) 



X5 = xe + 2xi 
x'g = -2x5 - 2x6 + 0.007xi -I- 0.6x2 

Xj = Xg — 0.5X2 

Xg = -5x7 - 4xg - 0.05x7 + 2.5x1, 

where again mg = 2, mi = 1 and fXoo = 3.8. 

Our theoretical results show that system (|7.54p . as well as the replicators, admit the chaos through 
period-doubling cascade and obey the universal behavior of Feigenbaum. Figure [TBI shows the 2— dimensional 
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Fi gure 15l 3 — dimensional projections of the chaotic attractor of the result-system 117.51I I. (a) Projection on the X3 — xq — xg 
space, (b) Projection on the xi — X5 — xio space. Since we are not able to provide a picture of the whole chaotic attractor in 
the 12— dimensional space, to obtain better impression, the projections shown in (a) and (6) are presented. These pictures give 
us an idea about the fascinating chaotic attractor of the result-system, which can possibly be called as 12D Chua's attractor. 



projections on the Xi — X2,X3 — 2:4,3:5 — xq and X7 — Xs planes of the trajectory of the result-relay-system 
(fTM)) with initial data a;i(0) = 1.37, X2(0) = -0.05, xaCO) 0.05, 2:4(0) = -0.1,X5(0) = 1.09,2:6(0) = 
—0.81,2:7(0) ~ 0.08,2:8(0) ~ 0.21. The picture seen in Figure [TBI (a) is the attractor of the generator (|7.53p 
and accordingly Figure 1161 (^) ~ (d) represents the attractors of the first, second and the third replicator 
systems, respectively. It can be easily verified that all replicators used inside the system ()7.54|) satisfy 
condition {A7). The resemblance of the chaotic attractors of the generator and the replicators is a consequence 
of morphogenesis of chaos. 

Next, we make use of the same trajectory to illustrate the 3— dimensional projections on the 2:3 — X5 — 2:7 
and 2:4 — xq — Xs spaces in Figure [TTI The visualizations shown in Figure [TTI (a) and (6), inform us about the 
shape of the result attractor of system (|7.54p . 

Now, let us continue with the control of morphogenesis of chaos. In paper [44], it is indicated that to 
control chaos of system (|7.53p generated through the period-doubling cascade, unstable periodic orbits of the 
logistic map /i(s,/ioo) — /^oos(l — s) must be necessarily controlled, and the OGY control method [5^ 1100) is 
proposed to stabilize the periodic solutions of system (|7.53p by favour of the conjugate control system 



Xi = X2 

x'2 = -0.182:2 - 22:1 - 0.000042:5* + 0.02cos(27rt) + i/(i,ti,/xO 



(7.55) 



In system (|7.55|) . the sequence is given by the formula 

where k'^^-', j > 0, with k^^^ = to G [0,1] is the target periodic orbit of the logistic map /i(s, /^oo) to be 
controlled, and k^+i = h{Ki, ^oo)i * > 0, with kq = ti G [0, 1]. 

In this method, the sequence jli is supposed to vary in a small range around the value jioo, that is, ^ can 
take values in the range [/ioo — ^, /ioo + S\, where (5 is a given small positive number. When the trajectory is 
outside the neighborhood of the target periodic orbit, no parameter perturbation is applied and the system 
evolves at its nominal parameter value ^00 ■ That is, we set /ii = /ioo when — /ioo| > 5. 
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Fi gure 161 2— dimensional projections of the chaotic attractor of the result-system lj7.54| l. (a) Projection on the xi — X2 plane, 
(b) Projection on the X3 — X4 plane, (c) Projection on the X5 — xq plane, (d) Projection on the — xg plane. The picture in 
(a) shows the attractor of the prior chaos produced by the generator l|7.53| l. which is a relay-system, and in (6) — (d) the chaotic 
attractors of the replicator systems are observable. The illustrations in (b) — (d) repeated the structure of the attractor shown 
in (a), and the mimicry between these pictures is an indicator of the morphogenesis of chaos. 



To simulate the result, let us consider the system 

x'2 = -0.18a-2 - 2xi - 0.00004a;? + 0.02 cos(27r<) + i^{t, ti,fl^) 
x'^ = X4 — O.lxi 

x'^ = —10x3 — 6x4 — 0.03a;| -I- 4x2 ry\ 
X5 = Xg + 2xi 

x'q = -2x5 - 2x6 + 0.007x1 -I- 0.6x2 

x'y = X8 — 0.5X2 

Xg = -5x7 - 4x8 - 0.05xy 4- 2.5xi 

which is the control system conjugate to the rcsult-relay-system (|7.54p . where again mp = 2 and mi = 1. 

To simulate the control results, we make use of the values 5 = 0.19, ti = 0.5, to — 2.8/3.8 and the trajectory 
of system (jT^T)) with the initial data xi(0) = 1.37,X2(0) = -0.05, X3(0) = 0.05, X4(0) = -0.1,X5(0) = 
1.09, X6(0) = -0.81, X7(0) = 0.08, X8(0) = 0.21. Taking the value = 2.8/3.8 means that the control 
mechanism is applied around the fixed point of the logistic map, and consequently stabilizes the 2— periodic 
solutions of the generator and the existing replicators. The control is switched on at t = 25 and switched 
off at t = 125. The graphs of the coordinates X3,X5 and x^ are shown in Figure [181 It is observable that 
the 2— periodic solutions of the replicators and hence of the result-relay-system (|7.54p is stabilized. In other 
words, the morphogenesis of chaos is controlled and the results of Theorem 16.11 are validated also through 
simulations. One can see that after approximately 60 iterations when the control is switched off, the chaos 
becomes dominant again and irregular motion reappears. 

7.4 Quasiperiodicity through morphogenesis 

Now, let us indicate that if there are more than one generator system, morphogenesis will have some new 
forms such as periodicity gives birth to quasiperiodicity. 
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Figure 17: 3— dimensional projections of the chaotic attractor of the result-rclay-systcm 117. 541 1. (a) Projection on the 2:3 —xs—xy 
space, (b) Projection on the X4 — xg — xg space. The illustrations presented in (a) and (b) inspire us about the shape of the 
impressive chaotic attractor in the 8— dimensional space. 



In paper |101| . it is shown that the Duffing equation 

x" + 0.168x' - f (1 - = fism{t) (7.58) 

where fj, is a parameter, admits the chaos through period-doubhng cascade at the accumulation point /ioo = 
0.21 of the sequence {/ij} of parameter values, that is, for fi = ^00, the equation (j7.58p admits infinitely many 
periodic solutions with periods 2fc7r, for each natural number k. Making the change of variables t = 2-ks and 
x{t) = y{s), one can obtain the equation 

y" + 0.3367ry' - 2TT'^y (l - y'^) = An'^fi sin(27rt). (7.59) 

Likewise equation (|7.58p . it is clear that the equation (|7.59p also admits the chaos through period-doubling 
cascade and admits infinitely many periodic solutions with periods 1,2,4,... when one considers the equation 
with the parameter value fi = ^co- Using the new variables xi = x^X2 = x' ^x^ = t and .T4 ~ y^x^ = y', xg — s, 
respectively, one can convert the equations (|7.58p and (|7.59|) to the systems 

x[ = X2 

4 = -0.168x2 + ^{l-xf)+ Aioo sin(a;3) (7.60) 
4 = 1 

and 

2^4 — ^5 

Xg = -0.3367ra;5 + 27r^X4(l - xl) + An'^fioo sin(27rx6) (7.61) 
< = 1 

respectively. Now, we shall make use both of the systems (|7.60|) and (|7.61|) as generator systems to obtain 
a chaotic system with infinitely many quasiperiodic solutions. We mean that the two systems admit in- 
commensurate periods and consequently their infiuence on the replicator will be quasiperiodic and one can 
expect that replicator will expose infinitely many quasiperiodic solutions. Consider now the 9— dimensional 



48 




50 100 150 200 

t 



Fi gur6 18l OGY control method applied to the result-relay-system l|7.54|l . (a) The graph of the xs — coordinate, (b) The graph 
of the x^— coordinate, (c) The graph of the xy— coordinate. For stabilizing the 2— periodic solution of the rcsult-rclay-system 
l|7.54|l . the OGY control method is applied to the generator system l|7.53|l around the fixed point 2.8/3.8 of the logistic map. 
The value 5 = 0.19 is used, and the control starts at t = 25 and ends at t = 125. After approximately 60 iterations when the 
control is switched off, the chaos becomes dominant again and irregular motion reappears. The presented coordinates 2:3, 2:5 
and XT, all correspond to replicator systems and it is possible to obtain similar illustrations for the remaining ones, which are 
not just simulated here. 

result-system 

x'r^ = -0.168x2 + 0.5x1(1 - xj) + ^00 sin(a;3) + Ci{x2{t - 2tt) - X2{t)) 
4 = 1 

X^ = X5 

x'^ = -O.SSGttxs -f 27r2x4(l - xf) + A^:'^^loo sin(27ra;6) + C2{x5{t - 1) - x^it)) (7.62) 
4 = 1 

Xy = + Xl + Xii 

Xg = -3X7 - 2X8 - 0.008X7 + X2+XZ 
Xg = -Xg + X3 -I- X6, 

where the last three equations are of a replicator system. 

We make use of the Pyragas control method together with the values of Ci = 0.1, and C2 = 0.62, to 
stabilize the 2n and 1 periodic solutions of the generators systems (|7.60|) and (|7.6ip . respectively. For our 
purposes, let us simulate a solution of system (|7.62p with initial data xi(0) = 0.4, X2(0) = — 0.1,X3(0) = 
0,X4(0) = -0.15, X5(0) = 0.4,X6(0) = 0,X7(0) = 1.1,X8(0) = 2.5,X9(0) = 0.12. The simulation results are 
seen in Figure [TH The control mechanism starts at t = 40 and ends at t = 120. Up to t = 40, the chaos 
not only in the generator systems, but also in the replicator is observable. Between i = 40 and t — 120 the 
quasiperiodic solution of the replicator is stablized. After t = 120, chaos in the system (|7.62p is develops 
again. 

Possibly the obtained simulation result and previous theoretical discussions can give a support to the 
idea of quasiperiodical cascade for the appearance of chaos which can be considered as a development of the 
popular period-doubling route to chaos. 

In paper |102j . it has been mentioned that, in general, in the place of countable set of periodic solutions 
to form chaos, one can take an uncountable collection of Poisson stable motions which are dense in a quasi- 
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Figure 19: Pyragas control method applied to the result-system 117.6211 . (a) The graph of the 0:2— coordinate, (b) The graph 
of the 3:5— coordinate, (c) The graph of the xg— coordinate. The simulation result for the result-system I I7.62I I is provided such 
that in (a) and (fe) periodic solutions with incommensurate periods 1 and 27r are controlled by Pyragas method and in (c) 
quasiperiodic solution of the replicator system is pictured. For the coordinates xi, X4 and we have the similar results which 
are not illustrated here. 



minimal set. This can be also observed in Horseshoe attractor |33| . These emphasize that om' simulation of 
quasiperiodic solutions can be considered as another evidence for the theoretical results. 

7.5 Replicators with nonnegative eigenvalues 

We recall that in our theoretical discussions all eigenvalues of the real valued constant matrix A, used in 
system (|2.12p . are assumed to have negative real parts. Now, as open problems from the theoretical point 
of view, we shall discuss through simulations the problem of morphogenesis of chaos in the case when the 
matrix A possesses an eigenvalue with positive or zero real part. 

First, we are going to concentrate on the case of the existence of an eigenvalue with positive real part. 
Let us make use of the Lorenz system (|5.36p together with the coefRcients cr = 10, r = 28 and 5 = | as the 
generator, which is known to be chaotic [351 167j . and set up the 6— dimensional result-system 

x[ = -lOxi + 10a;2 

X2 ~ —X2 + 28x1 — 2:13:^3 

2^3 — X1X2 — 3X3 ,„ 
x'^ = —2X4 + Xi 
= — 3X5 -I- X2 
x'q = 4X6 - Xg -I- X3. 

It is crucial to note that the system (|7.63p is of the form of system (|2.10p + (|2.12p where the matrix A admits the 
eigenvalues —2, —3 and 4, such that one of them is positive. We take into account the solution of system (|7.63p 
with the initial data xi(0) = -12.7,X2(0) = -8.5,X3(0) = 36.5, X4(0) = -3.4,X5(0) = -3.2,X6(0) = 3.7 and 
visualize in Figure [50] the projections of the corresponding trajectory on the Xi — X2 — X3 space and on the 
X4 — X5 — Xg space with three different standpoints. It is seen that the replicator system admits the chaos 
and the morphogenesis procedure works for system (|7.63p . 
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Fi gurc 20l 3— dimensional projections of the chaotic attractor of the result-system (17.6311 ■ (a) Projection on the x\ — X2 — x^ 
space, (b) Projection on the X4 — X5 — xq space, (c)-(d) Projections on the X4 — X5 — xq space with different standpoints. In (a), 
the famous Lorenz attractor produced by the generator system I I5.36I I with coefficients a = 10, r = 28 and 6 = ^ is shown. In 
(6) — (d), as in usual way, the projection of the chaotic attractor of the result-system l|7.63|l . which can separately be considered 
as a chaotic attractor, is presented with throe different standpoints. Possibly one can call the attractor of the result-system as 
6-D Lorenz attractor. 



Next, we continue to our discussion with the case of the existence of an eigenvalue with a zero real part. 
This time we consider the chaotic Rdssler system [67l I103| described by 

x'l = -(X2 + X3) 

4 = xi 0.2x2 (7.64) 

4 = 0.2 + a;3(a;i - 5.7) 

as the generator and constitute the result-system 

x'l = -{X2 + X3) 

x'2 ^ xi + O.2.X2 
4 = 0.2 + X3(a;i - 5.7) 
x'^ = —Axa + xi 

^5 = -X^ + X2 

4 



(7.65) 



-0.2x1 + X3. 



In this case, one can consider system (|7.65p as in the form of (|2.10p + (|2.12p where the matrix A is a diagonal 
matrix with entries —4,-1,0 on the diagonal and admits the number as an eigenvalue. We simulate the 
solution of system (|7.65p with the initial data xi(0) = 4.6,a:;2(0) = — 3.3,2:3(0) — 0,^4(0) = 1,0:5(0) = 
—3.7, xq{0) = 0.8. The projections of the trajectory on the xi — X2 — X3 and x^ — X5 — xq spaces arc seen in 
Figure [2TJ The simulation results confirm that the replicator mimics the complex behavior of the generator 
system. 

These results of the simulations request more detailed investigation which concern not only the theoretical 
existence of chaos, but also its resistance and stability. 



8 Conclusion 

Morphogenesis of different types of chaos such as obtained by period-doubling cascade, Devaney's and Li- 
Yorkc chaos is recognized in the paper. The definitions of chaotic sets as well as the hyperbolic sets of functions 
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a 



b 




Fi gurc 21; 3 — dimensional projections of the chaotic attractor of the result-system 117.65I1 ■ (a) Projection on the xi — X2 — ^3 
space, (b) Projection on the X4 — — xe space. The picture in (a) indicates the famous Rossler attractor produced by the 
generator system 117.6411 . The similarity between the illustrations presented in (a) and (b) supports the morphogenesis of chaos. 
The attractor of the result-system 17.6511 can be possibly called as 6D Rossler attractor. 

are introduced and the self-replication of the chaos is proved rigorously. The morphogenesis mechanism is 
based on a generating chaotic element inserted in a network of systems. Morphogenesis of intermittency 
as well as the double-scroll Chua's attractor and quasipcriodicity are discussed. Moreover, control problem 
of the self-replicated chaos is realized. Some of the results are illustrated through relay system's dynamics. 
Appropriate simulations are presented using the indicated method successively. 

One can construct the morphogenesis mechanism by the formation of consecutive replications of chaos 
or replication of chaos from a core system. It is also possible to construct a result-system using these two 
mechanisms in a mixed style. 

In the present paper, we use the term "morphogenesis" in the meaning of "processes creating forms" 
where we accept the form not only as a type of chaos, but also accompanying concepts as the structure 
of the chaotic attractor, its fractal dimension, form of the bifurcation diagram, the spectra of Lyapunov 
exponents, inheritance of intermittency, etc. Replication of a known type of chaos in systems with arbitrary 
large dimension is a significant consequence of our paper. More precisely, by the method presented, we show 
that a known type of chaos, such as period-doubling cascade, Devaney's or Li-Yorke chaos, can be extended 
to systems with arbitrary large dimension. We provide a new mechanism through the insertion of chaos 
from one system to another where the first is not necessarily a simple one as an "embryo" . To be more 
precise, we provide replication of chaos between unidirectionally coupled systems such that finally to obtain 
a result-system admitting the same type of chaos. 

Our method can be useful to explain extension of irregular behavior like crisis in economic models |66j as 
well as reproduction of chaos in mechanical systems |63j and electric circuits [97j . 
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